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AERODYNAMIC HEATING OF CONICAL ENTRY VEHICLES AT SPEEDS 


IN EXCESS OF EARTH PARABOLIC SPEED 

By H. Julian Allen, Alvin Seiff , 
and Warren Winovich 


SUMMARY 


The aerodynamic heating characteristics during Earth 1 s atmosphere entry at 
speeds greater than Earth parabolic speed are calculated for vehicles of conical 
shape. Ablative heat shields are assumed for these bodies and both laminar and 
turbulent boundary layers are considered. It is shown that if conical shape can 
be maintained, an optimum cone angle will exist and a cone of this angle will be 
superior to the more usual blunt entry shapes at speed well in excess of para- 
bolic speed. The improvements result from the fact that radiative heat-transfer 
contributions from the shock layer are much reduced for the cones so that, 
although the convective contributions are increased, a net gain is realized. For 
optimum cones, the approximate analysis indicates that the convective contribu- 
tions constitute 85 to 90 percent of the total heating. Solutions to the problem 
of maintaining conical shape as ablation progresses are considered, and some 
experimental demonstrations of means for accomplishing this are presented. 


INTRODUCTION 


Interplanetary travel may require entry into planetary atmospheres at speeds 
well in excess of parabolic speed in order to shorten trip times (ref. l) . Fig- 
ure 1 shows trip time as a function of entry speed into Earth 1 s atmosphere for 
travel from Mars and Venus . There are limits which must be set for the maximum 
entry speed allowable because of both the loads entailed (ref. 2) and the aero- 
dynamic heating. In this report we shall be concerned only with the heating 
aspects . 

Up to the present time, atmosphere entry speeds have been sufficiently low 
that aerodynamic heating has been essentially a convective process. Osborne 
Reynolds long ago (ref. 3) showed that the molecular process by which a fric- 
tional force is exerted on an aerodynamic surface is directly related to the 
process by which heat can be convected to that surface. In consequence, it can 
be shown (ref. b) that the incremental quantity of energy in the form of heat 
convected to an entry vehicle, SH, is related to the increment in time, St, by 



SH ~ FV St 


( 1 ) 


the proportionality 1 


where F is the total frictional force exerted on the vehicle. On the other 
hand, the incremental change in kinetic energy in the same interval is given by 

SE = DV St (2) 

where D is the total drag force experienced by the vehicle. It follows that 

SH ~ | SE (3) 

Let us consider that the ratio of the friction force to the drag remains 
essentially constant. The total energy in the form of heat convected to the 
vehicle during entry is then proportional to the total kinetic energy change 
•which occurs during entry, and, for a vehicle which is not to he destroyed on 
landing at the planet’s surface, this total kinetic energy is simply the kinetic 
energy of the vehicle at entry to the atmosphere. To minimize the convective 
heating, then, one must choose a vehicle shape with the smallest ratio of fric- 
tional force to total drag force. Thus one chooses "blunt shapes for which pres- 
sure drag is high, in order that the drag coefficient may he made as large as 
possible. The amount of heat transferred hy convection depends upon whether the 
boundary layer is laminar or turbulent. At the usual Reynolds numbers character- 
istic of entry, it is advantageous to have laminar flow. Accordingly, one 
strives to maintain such a flow. Also most entry bodies employ ablative shields 
for heat protection. In part, ablative systems are used because the ablating 
vapors fend off the air and so reduce the shear in the boundary layer from what 
it would be in the absence of these vapors, and, hence, the convective heat 
transfer as well. Under optimum conditions of high pressure drag and low fric- 
tional force, the fraction of the total kinetic energy convected to the entry 
vehicle can be kept very low indeed - of the order of l/lO of 1 percent or less. 

When, now, one considers the higher entry speeds desired in the future, one 
finds it difficult to prevent the heat transfer for a blunt body from increasing 
rapidly with increase in speed because convection is no longer the sole important 
heating mechanism involved. In the new speed regime the air which enters the bow 
shock layer undergoes such high molecular excitations that it becomes a powerful 
source of radiative energy. This process is well described in the literature 
(e.g., refs. 5 to 10). At this point it is only necessary to note that the most 
important radiative contribution to aerodynamic heating varies with velocity by 
as much as the fifteenth power and almost directly as the density. Hence, 
although for entry at near -Earth satellite speed the radiative contribution for 
a blunt body is usually trivial, it tends to become overwhelming at speeds well 
in excess of Earth parabolic speed. One is led, therefore, to re-examine the 
effect of vehicle shape when radiative contributions to heating are important. 

^All symbols are defined in appendix A. 
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To this end, consider the conical body shown in figure 2. As an entry vehicle, 
it has the disadvantage of its drag coefficient being less than that for an 
essentially flat -faced body so that, other factors being equal, it must accept a 
larger convective heat load. On the other hand, the radiation per unit of shock 
layer volume depends upon excitation of the gas in the shock layer and so depends 
upon the component of velocity normal to the bow shock. This radiation varies 
as about the fifteenth power of the sine of the shock angle. Thus, at suffi- 
ciently high entry speeds, a reduction a-n shock angle, which increases convective 
heat input but greatly reduces the radiative heat input, can provide a net gain. 
Under these conditions it is to be expected that an optimum cone angle exists for 
any given set of entry conditions. It is the purpose of this paper to analyze 
entry bodies of conical shape to find these optimums. The metric system of units 
(kilogram-meter -second) are employed throughout this analysis. 


ANALYSIS 


The rate of energy input in the form of heat to an entry body may be written 
(ref. 11 ) 


^=|c hP v 3 a (4) 

A number of assumptions will be made to simplify the analysis. The first 
two are that aerodynamic lift is zero and that during the time the heating 
process is important, the deceleration is large compared to the acceleration of 
gravity; hence, the gravitational effect can be ignored. In this case the 
trajectory is essentially a straight line (ref. 4 ) so that the time rate of 
change of altitude, Y, is 


— = -V sin 7 
dt 


(5) 


where y is the flight -path angle as measured down from the local horizontal and 
is constant. The third assumption is that the air density in the atmosphere 
varies exponentially with altitude (see ref. 4, 11 , or 12 ). 


-PY 

p = pp 0 = p Q e 


(6) 


where p is the air density in terms of a reference sea- level density, p Q , and 
P is a constant. 

From equations (5) and ( 6 ) then 


dt 


dY = dp 
V sin / ppV sin 7 


(7) 
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The fourth and fifth assumptions are that the drag coefficient is constant, 
and that the mass remains essentially constant during the entry. Thus, from 
reference k, 


V 2 = V E 2 e" Bp (8) 

where V E is the vehicle speed at atmosphere entry and B is the ballistic 
coefficient 


B = 


C D P 0 A 
pm sin 7 


(9) 


■where m is the entry hody mass. 

With equations (4), ( 7 ), and (8), the energy input to the entry vehicle in 
the form of heat for the whole entry is 


H = 


Pq V E A 
2 P sin 7 Jp =0 


C H e" Bp dp 


( 10 ) 


The total kinetic energy change during entry is 


E = | mV/ 


( 11 ) 


since the final speed at landing is zero for a vehicle which is to land intact. 
Then equations ( 9 ), (10), and (ll) combined give the fraction of the total 
kinetic energy which appears as heat to the vehicle, tj, as 

H = 5 .JL r o H e - B5 dp (12) 

® J p=0 H 

The sixth assumption is that the heat shield is the ablation type for which 
mass loss is only by the process of sublimation or vaporization. Let £ be the 
heat required to bring a unit mass of ablator from the cold state through vapori- 
zation expressed in kinetic energy units (i.e., in units of square of velocity). 
Then the mass loss by ablation in terms of the entry vehicle mass will be 


£ - » ® <13) 

We now proceed to determine the energy fractions, tj, as a function of entry speed 
for conical bodies under the assumption (seventh) that the cone angle of the con- 
ical body remains unchanged during ablation. To this end we determine tj as the 
sum of the contributions due to equilibrium radiative heating (denoted by rj e ), 
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the none qu il lb r ium radiative heating (denoted hy T] n ), and either the laminar 
convective heating (denoted hy tj^) or the turbulent convective heating (denoted 
hy r]^). The eighth assumption is that each contribution may he calculated 
independently of the other. Cases for mixed laminar flow and turbulent flow are 
not treated hut the limiting effect of approach to free -molecular flow is 
included. The drag coefficient used in the evaluation of tj is assumed to he 
the Newtonian value (ninth assumption) 


C D = 2 sin 2 Q C 


(u) 


where 0 q is the half -cone angle. Thus, base pressure is ignored since the 
speeds of interest are great, and effects of friction on drag are ignored on the 
presumption that the cones of interest are never slender enough to warrant the 
complication of including this effect. 


Equilibrium Radiative Heating 


For the purposes of this analysis it has been assumed (tenth and eleventh 
assumptions) that the radiation per unit volume from a shock layer which is in 
thermodynamic and chemical equilibrium is constant throughout the shock layer 
(i.e., energy depletion due to radiation is ignored) and is nonabsorptive within 
the layer. The equilibrium radiation (appendix B) is, under these conditions, 
determined solely by the density, p, ahead of the bow shock and the velocity 
normal to it. 


u = V sin e w (15) 

An examination of the available data indicates (see appendix B) that the 
time rate of equilibrium radiation per unit volume of gas cap can be expressed 
approximately for the lower speed range 


dEe 

dv 


C 



1 


p p 


U < Ux ^2 = 13,700 m/sec 


and for the upper speed range 


(l6a) 


dE e 

dv 


= C e£ U 





U > Ux j 2 = 1-3,700 m/sec 


(16b) 


wherein 
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p = 


'll = 


q . 2 

Ce x 

C e 2 


1.80 

15.45 

5.05 


= 6.i4xl0- 49 kg 

m sec 3 4 m 


15 . 45 


(IT) 


= 6.WBOO-' -^6-5 

m sec 3 V m 


5.05 


The time rate of the total equilibrium radiation from the shock layer is 


E e = 


cffi e 

dv 


(18) 


and approximately half of this radiation is received hy the body surface , if 
there is no surface reflectivity. Then the equilibrium-radiation heat -transfer 
coefficient is 


_ (l/2)E e _ (<lEe/dv)v 
CH e " ( l/2 ) P V 3 A ” p Q pV 3 A 

where 

/tan 2 0 W - tan 2 0 C \ 

V = Jtr ' b \ 3 tan 3 0 q ) ^ 

and 

A = rtr^ 2 

J 

while for the ARDC standard atmosphere (ref. 13) 

p Q = 1.225 kg/m 3 

so that from equations (l6), (8), and (14), equation (12) becomes 2 for 
Vj; < ©W 


(19) 


( 20 ) 


( 21 ) 


2 Properly, the value of C He should be made zero in the free -molecular 

flow regime but the strong dependence of T] e on p in equations (22a) and (22b) 
mak es the contribution of the integral to Tj e for the free-molecular altitudes 
negligible for vehicle weights and sizes of usual interest. 
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I 


<3.1 -3 

Tig = Br b c ei v e 


sin 1 %(tan 2 % - tan 2 0 q) 


6 sin 2 0 Q tan 3 0 q 
and for Vg > U lfS /sin. % 


T 1 -(B^K^-Op^^ 

/ e dp 

; p=o 


( 22 a) 


'll " 3 

T!e = Br^ceje 


sin 1 0 ^(tan 2 % - tan 2 e c )l n 1 -(B/ 2 )(q -i)p -n., 

/ e P dp 


+ Bi^ce^e 


q.p -3 


6 sin 2 00 tan 3 00 
sin ^ 2 0 w (tan 2 0 W - tan 2 0 r ) 


P=P 


1,2 


6 sin 0 q tan 3 0 q 


' P p l ,2 -(B/ 2 )(q -i)p 

e p p x dp 

(22b) 


'P=o 


where 


p = — Zn 

^ 2 B V U x , 


E sin 0 W 


c e , = ^ = 5-012X10" 49 ( — 


sec 3 V m 


ce 2 - ^ - 5.318X10-= JSL 
2 P 0 sec 3 \ m 


15 • 45 


5 • 05 


(23) 


Details of the integrations of equations (22a) and (22b) are given in 
appendix B. The solution may be given in the form 

r b*e Z e 


^e B°‘ s 


(24) 


wherein is a function of 0 q (see appendix B and table I) and Z e is a 

function of V e sin 0^ (see appendix B and table II). 


Nonequilibrium Radiative Heating 


For the purposes of this analysis it has been supposed that the portion of 
the shock layer which is not in thermodynamic and chemical equilibrium is con- 
fined to a region so close to the bow shock wave that it may be considered to 
originate at and be proportional to the area of the bow shock wave. In addition, 
the nonequilibrium radiation process is regarded as one involving binary colli- 
sions so that the radiation is independent of air density until "collision 
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limiting" occurs (see ref. 10 and appendix C) and is a function only of the 
velocity normal to the shock -wave (see eq. (15)). Collision limiting is presumed 
to begin when p = p C £, in th e regime for which p < p c ^ the nonequilibrium 

radiation is considered to vary directly with density. When the free -molecular 
flow regime is reached (i.e., when p < ), this radiation is set to be zero. 

Thus, the time rate of nonequilibrium radiation per unit area of bow shock wave 
is given as 


dEn q 

aT = c * u ; 5 > p cZ 


dEn 
da P, 




cl 


pU S ; P C Z > P > Pfm 


dE- 


'n 


da * ° 5 5 < 5 fm 

An examination of the available data indicated (see appendix C) that 

C n = 0.7 4x10 " 22 sec 4 /m 4 

s = 7 

PcZ = 10-3 


( 25 a) 


( 25 b) 

( 25 c) 


(26) 


The total time rate of radiation from the gas cap which is not in equilibrium is 
then 


E 


n 


dE n 

a 

da 


(27) 


and only half of this radiation is received at the vehicle surface if the surface 
reflectivity is zero. Then the nonequilibrium radiative heat -transfer coeffi- 
cient is 


C H n 


(l/2)E n _ (dE n /da)a 

( l/2 ) pV 3 A ' p o pV 3 A 


(28) 
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•where 


a = otiv 


A = rtV 


tan 2 


in % tan^ 0 q. 


while, again. 


p n = 1.225 W^ 3 


(29) 


so that from equations ( 25 ), (8), and (l4), equation (12) hecomes 


s -3 /si n 5 ' 1 6 W tan 2 

\2 sin 2 0 q tan 2 0 q. 


1 n = Bc n V E 


^ -(B/2)(b-i)p 

e dp 


J P=P 


cl 


+ - 
P, 


' P cl -(B/2)(s-l)p 


ap 


c Z a p=p 


fm 


(30) 


where 


; n = — = 0.60X10 " 22 sec 4 /m 4 


(31) 


Details of the integration of equation ( 30 ) are given in appendix C. The 
solution may he given in the form 


/v E y 
= Y n®n ( j 


(32) 


where is a function of 0 q (see appendix C and table I) and 0 n is a 

function of B (see appendix C and table III). 


Laminar Convective Heat Transfer 


Calculation of the laminar convective heat -transfer coefficients as a 
function of velocity, air density, and cone angle was made by the procedure 
given in appendix D. In the absence of ablation it is indicated that the heat- 
transfer coefficient may be expressed as (sub zero indicates no ablation) 


% 


o 


C Z 


(33) 
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where c^ and j are functions of the -velocity V. When vapor ablation occurs, 
the vapor layer fends off the air and so reduces the convective heating. It is 
assumed in this analysis that for ablation (twelfth assumption) 


% 


1 + K^V 2 




(3*0 


where a ^ is an asymptotic lower limit when 

K z V 2 -*°° (35) 

and Ki is a constant depending upon the ablative material and can be defined as 


Oj-J 

K Z = j (36) 

wherein a z depends upon the molecular weight of the ablating vapors 3 and £ is 
the heat energy per unit mass (expressed in kinetic energy units of square of 
velocity) required to heat the ablator from the cold state through vaporization. 

As an approximation for speeds up to 13,000 m/sec, the variable c z and 
the j values may be replaced by the constants (see appendix D) 


c z = c Zl = 3. 5x10 ~ 4 sin 20 q , m 1//2 

t " 0 


(37) 


and at speeds above 26,000 m/sec, the variables 
constants 

c 2 = c l £ = 19*8 s] sin 20 q , 
j 2 = 1.17 

The foregoing applies in a continuum flow regime. In the free -molecular 
flow regime, the fraction of the total kinetic energy which, converted to heat, 
appears as heat to the vehicle is assumed to be one-half (ref. 15)j that is to 
say, in the free -molecular flow regime 


can be approximated by the 


m 


1/2 f ^ 


sec 


1. IT 


(38) 


% 



( 39 ) 


3 The value of 0,2 is l/2 3 l as defined in reference 14. 


10 



Of course the continuum results (eqs. (33) and (34)) only apply when the 
Knudsen number (ratio of air mean free path to body diameter) is very small com- 
pared to unity; conversely, the free -molecular result (eq. (39)) only applies 
when the Knudsen number is very large compared to unity. A transition from one 
regime to the other occurs when the Knudsen number is of order unity. Experiment 
has indicated that the transition from one regime to the other is smooth (see, 
e.g., ref. 1 6). For the present calculations it was deemed adequate to consider 
that the continuum result applied at all air densities for which 



( 1 - , ' 
VI + k 7 v 2 + a \ 


< 


- sm* 


0C 


(1+0) 


Within the free -molecular flow range, for the sizes and masses of entry vehicles 
of interest in this analysis, the velocity can safely be taken to be the entry 
velocity, V E . Thus the continuum flow results will be applied for all air den- 
sity ratios greater than 


fm 




LVi + I- 


+ a- 


V E J -vAb sin2 e C~ 


(41) 


with 


il = K Z V E 2 

and in the free -molecular flow range we take 


C D 

c Hz = — = sin 2 e c ; 


P < P. 


fm 


Finally, then, the laminar convective energy fraction may be written 
Bc z (l - a z ) p 1 e -B[x-(j/ 2 )]p 


^ l = 


ap 


2V E J ^ sin 2 e c J p fm (1 + | ie _B P) 

Bc l a l 


2V E j N /r^ sin 2 9 C A 


1 e -B[i-( J/s)]p a.p + B 

2 > 


fm 


fm 




e' B ^ dp 


(42) 


(43) 


(44) 


It is shown in appendix D that the solution of equation (44) for the entry 
speed range up to 13,000 m/sec may be given as 
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( 45 ) 


’U - ^Zx = *hjk c(l - a Z^Zx + 1 -772a z - ® Zl ] - 2sr z 


B 


" a l 


Fb/ VI + £ 


+ a- 


and, for the entry speed range above 26,000 m/sec, as 


¥ 


^ = ^Z 2 = 


Z£ 


(V E /l0 4 ) 


4sl.1T j r - 


~ [(1 


" °l^i z + 2 *748a^ 


- v 


- 2 


*Z; 


(v E /io 4 r 


B 


4s£.S4 V.1 + £ 


- a-, 




(46) 


where 

functions of e c (see appendix D and table I) 

functions of (see appendix D and table IV) 

functions of B (see appendix D and table III) 

For speeds between 13,000 and 26,000 m/sec the values of are assumed to 

change smoothly from the first range to the second, and the following applies: 

Define 


Ox = r|j 1 at V E = 13,000 m/sec 
G 2 = T] 2 2 at V E = 26,000 m/sec 

G3 = (t)j at V E = 13,000 m/sec) - (tj z at V E = 12,000 m/sec) 
G4 = (n at V E = 27,000 m/sec) - (r ) z at V E = 26,000 m/sec) 
Then calculate 

b 0 = Ox + 5(G 2 - Gx) - 52G 3 - 26g 4 


} (47) 


*>i = - ( yj )(G a - Gx) + 8OG3 + 5OG4 


*2 = (?f?)( g £ - G i) - (^r) g 3 - ( g 4 


b 3 = - 


169/ 

2000 

2197. 


13 


400 

13 




(a. - on + ( m) 03 ♦ '***> 


169 


VI69 


g 4 


( 48 ) 
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so that for the speed range 13,000 m/sec < Vg < 26,000 m/sec 


n* = ^i x . £ = ^o + 



( 49 ) 


Turbulent Convective Heat Transfer 


The analysis of appendix E shows that in the absence of ablation the 
turbulent heat-transfer coefficient can be expressed with reasonable accuracy by 



9q cos 

- 

sin 6 ^ sin 6 q 


c t 


■V _0 . 204 0.148 

i y. 




for 0 q expressed in degrees of arc and, 
for V sin 0^ 6 7,500 m/sec 

c t = c-t 1 = 8.80xl0“ 7 m°* 148 (deg) _1 (m/sec) _ °' 66 1 

k = ki = -0.66 j 

and for V sin 0w = 7,500 m/sec 

c t = c t 2 = 1 -35xl0" 3 m°* 148 (deg) -1 (m/sec) 0, 16 
k = k 2 = 0.16 


(50) 


(51a) 


(51b) 


The effect of vapor ablation to reduce the turbulent convective heat transfer can 
be assumed in the form (twelfth assumption) 


CH t = U + °t) ^ 


K - ttt 

K t -T 


(52) 


wherein depends upon the molecular weight of the ablating vapors." 


4 The value of ext = l/2 Pt as defined in reference 14. As suggested in the 
reference, we shall use ext =* l/3 in this report. 
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As in the case of laminar flow, one properly should consider separately the 
free -molecular flow regime. However, the interesting turbulent flow cases are 
those for which B is small (no more than a few hundred at most). With this 
restriction the presumption that the continuum regime extends up to entry con- 
ditions (p = 0) in calculating the heat transfer leads to negligible error. 

Then equation (12) gives for 


V E = " I — ° W sec 


sin & 


W 


= Bc t! 

Tlt vJ * 1 * 0 - 148 


00 COS 00 


. ki 


1 - G b 


E r b ^ sin % sin 0 C 7 u p=o XL + 


-Bp 


+ at 


- b G4> 


dp 


-0.204 


(53a) 


while for 


, r > 7500 , 

Vtt - —— m/sec 

* “ sin ft 1 


W 


Bc ti 


0Q COS 00 


V-g Kl r^° * 148 ^2 Sin kl " '"•- 3 


1 - G t 


1 % slna 2 , T / V E sin V l+ l + e _B P V p°-2°4 

^ B 7500 ') 


e -B(~¥> 4 . 

e dp 


Bc + 


k 2„ 0.148 \o ^ k 2 


0Q COS 0Q 


/v E sm^ 

, B \ 7500 J 


V E “^r b u * J -‘ i0 ^2 sin“ i; 0 w sin 3 0(/^p=o 


1 - 


1+ £ + e 


g t , g V V ^ dp 

■Bp U t J .0.204 


(53 B ) 


wherein 


St = W 


(54) 


It is shown in appendix E that the energy fraction may he given in the form 


B t = 


B 


0.204 


O. 148 


At " ®t*t 



(55) 


wherein is a function of 00 (see appendix E and table I), is a function 

of B (see appendix E and table III), and, for given values of a b and Kt? 
a function of both Vg and 9q (see appendix E and table V). The second term in 
equation (55) will only be important when the ballistic parameter,' B, is very 
small. 
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RESULTS 


The analysis of the preceding section was used to compute by IBM 7090 the 
heat-transfer characteristics of conical entry vehicles entering the Earth 1 s 
atmosphere. Since the analysis requires attached bow shock waves which limit 
the maximum cone half -angle to about 55° * results at larger angles are shown by 
dotted curves in all the figures to follow. The calculations have been carried 
out for a 1-meter base radius. In one case the heat shield was assumed to be 
composed of a low temperature ablator having the assumed characteristics of sub- 
liming Teflon, and a high temperature ablator having the assumed characteristics 
of vaporizing quartz. For the Teflon the characteristics assumed are 


t, = 2.2x10 s m 2 /sec 2 


Kj = 12xl0“ 8 sec 2 /m 2 (for 

= 0 . 26 ) 

KL(. = 4x10“ 3 sec 2 /m 2 (for 

0 ^ = 0.09) 

a l = <*t = °- 1 



(56) 


For the quartz the characteristics assumed are 

£ = l6xl0 6 m 2 /sec 2 

= 1.5xlG~ 8 sec 2 /m 2 (for = 0.24) 

> (57) 

K-j-. = 0 . 5x10"“ 8 sec 2 /m 2 (for cx^ = 0.08) 
a l = °t = °> 1 

J 

The choices of a a-^, and £ were based upon available literature (e.g., 
refs. 14 and 17 ) . The choices for cr^ and crt are based on the knowledge from 
experiments with Teflon (ref. 18) that the asymptote is not zero (as predicted 
by the usual theories). The asymptotic value is probably a function of molecular 
weight and thermal conductivity at least, and is probably different for laminar 
and turbulent flow. Values from experiment range from 0.05 to 0.2. The choice 
of 0.1 is an arbitrary one. 


Laminar Flow 


The analysis for all laminar flow during vehicle descent (excepting free- 
molecule flow at entry) is probably restricted to those cases for which the max- 
imum Reynolds number does not reach too large values. This condition is only 
fulfilled for the larger values of the ballistic parameter (B of the order of 
hundreds or thousands) and the results presented are so restricted. It is well 
to treat, first, a single case to illustrate the typical effect of increasing 
entry speed on the variation of the contribution of radiative and convective 
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heating to the total. Such an example is given in figure 3 tor which the entry 
hody has a cone half -angle of 30°, a “ballistic parameter of 200, and a Teflon 
head shield. At low entry speeds the convective transfer dominates and in this 
range the total energy fraction falls with increase in entry speed because 
(eg.. (34)) the effectiveness of ablation to reduce the convective heat-transfer 
rate is assumed to improve with increasing speed and because Ch^ 0 diminishes 
(above 10 km/sec) with increasing speed. At the higher entry speeds, the radia- 
tive heating rises rapidly with speed so that the total energy fraction exhibits 
a minimum. Note that the contribution of the nonequilibrium radiation to the 
total energy fraction is very small and could well be ignored. This observation 
applies generally to all cases of interest. 

Consider next the more general case wherein the ballistic parameter is fixed 
but we vary the cone half -angle over a wide range. The energy fractions as a 
function of entry speed are then typically those shown in figure 4 for, again, a 
ballistic parameter of 200 and a Teflon heat shield. The envelope values giving 
the minimum energy fraction as a function of entry speed are shown by the dashed 
curve . 

If now we consider various values of B appropriate to laminar flow, we 
can determine a series of envelope curves. Such envelopes are shown in fig- 
ure 5 (u) for a Teflon heat shield and in figure 5(b) for a vaporizing quartz 
heat shield. It is seen that, at any entry speed, lowering the ballistic coeffi- 
cient diminishes the energy fraction. 


Turbulent Flow 


The analysis for the turbulent flow case made no allowance for any laminar 
flow which is, of course, unrealistic. However, It is expected that the errors 
resulting from failure to allow for any laminar flow will be unimportant if the 
maximum Reynolds number occurring during the flight trajectory is very high. 

Thus we can presume that the "all turbulent" results will be realistic if the 
ballistic parameter is small - say of the order of 20 or less - but will be In 
error by unknown magnitudes for larger values. Accordingly, the results pre- 
sented are restricted to the smaller ballistic parameters. As for the laminar 
case it is well, first, to illustrate the typical effect of increasing entry 
speed on the variations of the contribution of radiative and convective heating 
to the total. Such an example is given in figure 6 for which the entry body has 
a cone half -angle of 30°, a ballistic parameter of 20, and a Teflon heat shield. 
The variation of the energy function with entry speed is seen to be similar to 
that for the laminar flow case (fig. 3) and for the same reasons. Again, it 
should be noted that the contribution of nonequilibrium radiation to the total 
is very small and could well be ignored. As for the laminar case, this observa- 
tion applies generally to all cases of interest. 

Consider, next, the more general case wherein the ballistic parameter is 
fixed but we vary the cone half -angle over a wide range. The energy fractions 
as a function of entry speed are typically those shown in figure 7 for a 
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■ballistic parameter of 20 and a Teflon heat shield. As "before, the envelope 
values are given hy the dashed curve. 

If, now, we consider various values of B which may he appropriate to 
turbulent flow, we can determine a series of envelope curves. Such envelopes 
are shown in figure 8(a) for a Teflon heat shield and in figure 8(b) for a 
vaporizing quartz heat shield. It is seen that at any entry speed, the energy 
fraction is again diminished by lowering the ballistic coefficient. 


DISCUSSION 


In the following, the salient features of the results for the laminar and 
turbulent flow cases are considered separately and then compared. Finally, the 
important assumptions made in the analysis are reviewed to assess their adequacy. 


Laminar Flow 


In figures 5(a) and 5(t>) it is clear that the least energy fraction is 
obtained by making the ballistic parameter as small as possible. Since the maxi- 
mum flight Reynolds number varies inversely with the ballistic parameter (see 
appendix F) and since there is probably an upper limit to the Reynolds number 
which one can allow and still enjoy laminar flow for the whole of the trajectory, 
then it follows that there is some minimum ballistic parameter at any given entry 
speed which can be permitted if the assumption of laminar flow over the whole 
trajectory is to apply. 

At lower supersonic speeds than we consider here, there is some support for 
the contention that there is some limiting Reynolds number (denoted hereinafter 
as ReQ_j_ m ) above which one cannot expect to maintain laminar flow. It is 
probable that in this high-speed region varies with free -stream enthalpy 

and with the composition of the ablation material. 5 

We have assumed from experience with nonablating surfaces at lower speeds 
that a maximum Reynolds number based on local surface flow conditions of 10 7 can 
be reached before turbulence occurs. However, we will show the effect of increas- 
ing or decreasing this limiting Reynolds number by presenting results for limit 
Reynolds numbers of 2xl0 7 and 0.5X10 7 as well. The analysis of appendix F gives 
the relation between the cone half -angle, the ballistic parameter, and the entry 
S p ee & for an arbitrary limit Reynolds number. The optimum values of ballistic 

^Not only is it probable that this limiting Reynolds number will be influ- 
enced by the molecular weight of the ablated vapor but it is most likely to 
depend upon the uniformity of the ablative process. For example, one expects a 
composite ablator, such as a plastic impregnated fibrous material or a charring 
ablator, by virtue of the jet-like injection of vapor from such a surface, to 
behave differently than some ablator such as uniform Teflon which should be free 
of such jetting. 
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parameter, cone half-angle, and energy fraction as a function of entry speed were 
found in the following way: For a series of entry speeds, Vg, "base radii, r^, 

ablation asymptotes, Oi , and ablation constants, , energy fraction, rj, was 
plotted as a function of cone half -angle, Qq , for various values of the ballistic 
parameter, B. (An example is shown by the solid curves in fig. 9*) On these same 
plots the values of B and 0c corresponding to limit Reynolds numbers of 0 . 5x10 7 , 
10 7 , and 2xl0 7 were located and the corresponding curves of r\ as a function of 
0 q were constructed from which one can determine the optimum values of energy 
fraction, il 0 p^ and the corresponding values of ballistic parameter, Boptj and. 

cone angle, 0c O pt‘ (For the particular case given in figure 9, these are shown 
by the dotted curves. The lowest values of t] correspond to r^pp . ) From the 
total complex of plots, it was then determined that for the range of variables 


0.5X10 7 < Re lim 2x10 7 

0.2111=^ = 1111 

0.05 = ^ 0.2 


2xl0 " 8 i§r = K z = i5xio -8 


sec 2 

m2 








the optimum could he expressed approximately hy 


V - Vt( Re limXl 0 - 7 )' O ' 281 r b O ' O44 (l0a z ) O * 31O (K z xl0B) 


-0*276 


«Copt = 5 0 opt (Ba lim XlO -V)-°-°8S rb 


Y ^ -O . 08 Siv .o . 0 17 ( 10(Jj ) O . 020 (K!X10 3 ) -O . 023 


-) * (~) 

r b/opt \ r b/ opt 


/ y . • <_><_> X r\A-\f \ ”0 • 033 > Q, Oi 035 

( Re lim x10 ” 7 ) r b °* O 41 (l 0 a 2 ) (K-jXlO 8 ) 


(58) 


(59) 

(60) 

(61) 


where the values of r] 0 pp, ©Copt’ an< l (®/ r b)opt are Fhe Functions of entry speed 
given in figure 10 and correspond by definition to values of r^p, Qq t , and 
( B /r-p ) opt when P p 


® e lim ~ P( ^ 7 
r b = 1 m 

a l = °- 1 

= 10 -8 sec 2 /m 2 




(62) 
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A fact of considerable importance is that when the energy fraction is 
optimum, the convective contribution is the principal part of the total, approxi- 
mately O .87 ± 0 . 02 . Thus it is important to know the convective contributions 
with an order of magnitude greater accuracy than the equilibrium radiation. The 
nonequilibrium radiation is, as noted earlier, a trivial contribution to the 
total. Other points to note about these optima are the following: 

(l) The optimum energy fraction is insensitive to the vehicle size (r^,) and 
can be approximately given as 


~ - ( a i \°' 2 3 - ( a it \°’ 3 

lopt = V 'vKjBeuJ ■ lopt W<W 


( 63 ) 


so that other factors being equal the optimum energy fraction is unfavorably 
affected by improvement in the heat of ablation as shown in figure 11 , as would 
be expected (eqs. ( 3 ^) and ( 36 )). 

The optimum mass loss ratio is given (eq. (13)) by the approximation 



^opt 


2 £ 


^ opt^E 


a Z Re lii 


O. 3 


-0.7 


( 64 ) 


The parameter 




(65) 


is plotted in figure 12 as a function of Vjj. 

From the foregoing it is seen that although vaporizing quartz, a high- 
temperature ablator, has a heat of ablation about eight times that for Teflon, 
a low-temperature ablator, the mass loss for quartz is only about one -fourth 
that for Teflon. These materials are compared in figure 13 wherein a limit 
Reynolds number of 10 Y is assumed. Note that at the optimum conditions one can 
keep the mass loss to the order of 10 percent or less for speeds well in excess 
of escape speed. 

(2) With reference to the optimum cone half -angle, we note that it is 
essentially independent of all factors except the entry speed and that 


^opt ~ ^opt 


( 66 ) 


As seen in figure 10, then, the optimum conical bodies are not slender ones. 
Even at an entry speed of 30 km/sec, the optimum cone half -angle is 25° of arc. 
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(3) Concerning the optimum ballistic parameter we note that it varies almost 
directly with the base radius, nearly inversely with the limiting Reynolds number, 
and is essentially independent of other factors except entry speed. Thus 


(67) 


The fact that B/r^ appears as the basic variable should not be surprising. As 
noted earlier nonequilibrium radiation is a trivial contributor to the total 
heating. Equilibrium radiation (eq. (24)) variation is 


B 


' e B °- 8 V r b. 


-o. a 


(68) 


while laminar convection (eq. (45) or (46)) is directly determined by B/r-^. 

With reference to the (B/r^) th 
made regarding trajectory angles. Since 


With reference to the (B/z^) there are some important connotations to be 


and since 


B 


opt ~ 


C P p o A 
Pm sin y J opt 


2p Q A sin 2 0 Copt 
Pm sin 7 


(69) 


while 


^opt ^opt 


m = 


rtr b 3p b 


3 tan 0 c opt 

where is the average density of the entry body, then for 


p Q = 1.225 kg/m 3 

b = — m -i 


7000 


we get 




r b/ 


= 5. 2x10 4 


opt 


sin2 ®C op t tan 
sin 7 


opt 


(70) 


(71) 


(72) 


(73) 
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or, at optimum. 


PfcT-h 2 sin 7 S 5.2x10“* 


sin 0 C O pt tan e Copt' 

WAbW 


(E^Xl 0 - 7 ) 


0.9 


(74) 


and the right side of this expression is a function only of the entry velocity 
and Reynolds number limit. The relation is shown in figure 14. Let us now 
review values of the product p^r.^ 2 sin 7 which are typical of entry vehicles. 

For current manned vehicles the order of values for density, base radius, 
and trajectory angle is 

= 250 kg/m 3 
r b = 1 m 
sin 7 = 0.06 

so that 

P^ 2 sin 7 = 15 kg/m (76) 



so that one expects that over tjae likely range of entry speeds permitted by entry 
load considerations - say, less than 20 km/sec (see ref. 2) - there may not be 
difficulty in maintaining laminar flow in the boundary layers. Of course, the 
small values of 7 which are required from load considerations demand the use of 
lift in order to hold the vehicle in an earth curvature flight path. Half -cone 
bodies (see, e.g., ref. 19) suggest themselves for such applications. Moreover, 
it should be noted that the analysis given herein for the aerodynamic heating of 
conical bodies is predicated upon the assumption of constant 7 trajectories 
and, hence, can only be considered a crude approximation for shallow-angle 
entries . 

For unmanned space probes the vehicle density will be about that for manned 
vehicles but the size will generally be less (say, r^ = 0.5 m), and the vehicle 
load consideration will usually permit steep descent (i.e., sin 7 up to unity). 
Steep descent removes the necessity for employing lift during entry while freedom 
in the choice of 7 may greatly ease the guidance problem prior to entry. How- 
ever, if we intend to maintain laminar flow over such vehicles, we may be forced 
to resort to flat trajectories except for low entry speeds, as seen by the 
results shown on figure 15 which give the maximum entry angle as a function of 
entry speed for typical vehicle quantities (pj., = 250 kg/m 3 , r b = 0.5 m) . As 
vehicle size increases, the probable trend with time, the choice of entry angle 
can become even more stringent. Thus, even for instrument probes, we may desire 
flat trajectories at the higher speeds - but not because of load restrictions. 

It must be remembered, however, that these conclusions are based upon the dubious 
assumption that the maximum Reynolds number permitted for laminar flow is of the 
order of 10 7 . The need for high-speed laminar -flow research is clearly evident. 
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Turbulent Flow 


As for laminar flow if one is to minimize the energy fraction one must 
employ the lowest ballistic parameter that can be permitted, particularly at the 
higher speeds (see figs. 8(a) and 8(b)). However, it should be noted at the out- 
set that one cannot indefinitely reduce the ballistic parameter and still recover 
the vehicle intact. In the absence of auxiliary drag devices (drag brakes, para- 
chutes, etc.) for terminal deceleration, the impadt speed (see eq.. (8)) is 

V 0 = V E e" B/2 (TT) 

so that values of B upward of 10 mast he set as a lower limit. With auxiliary 
drag devices, this lower limit might perhaps he halved. Even at these very low 
values of the ballistic parameter the turbulent energy fractions (see figs. 8(a) 
and 8(h)) are not as low as can he attained with laminar flow for limit Reynolds 
numbers of the order of 10 (l) * * * * * 7 . For the turbulent case, as for the laminar, when 
conditions are optimum, the convective heating is the dominant contribution to 
the energy fraction. Again, the optimum cone half -angles are not small even at 
very high speeds. 

In comparing the optimum energy fractions for turbulent flow with laminar 
flow it should be noted that the turbulent values are an order of magnitude 
higher (the ordinates of figs. 8(a) and 8(b) are expressed in percent while those 
of figs, 5(a) and 5(b) are tenths of percent). This disparity is forcefully 
demonstrated in figure 1 6 wherein for three entry speeds, optimum energy fraction 
is plotted as a function of the ballistic parameter. (Note the lower ends of the 
laminar curves represent a Reynolds number limit of 10 7 and that the dotted 
curves are arbitrary fairings to indicate how the turbulent and laminar might 
join. ) It is clear that laminar flow is to be sought even if the limit Reynolds 
numbers are low. Figure 17 shows the corresponding variations of optimum cone 
half -angles with entry speed and ballistic parameter and again emphasizes that 
the optimum angle is only, in essence, dependent on velocity. 


Review of Assumptions 


In the analytical development many assumptions were made. While a number of 
these are clearly valid approximations there are many that should be reviewed to 
ascertain, in retrospect, their adequacy. Discussion of them is given in the 
following : 

(l) The assumption that flight path angle is constant, which was made for 

the purpose of simplifying the analysis, is generally admissible when the entry 

trajectory is steep (see ref. 12). However, for manned vehicles and for high- 

speed probes which may enjoy laminar flow the assumption is not strictly admis- 

sible so that this restriction must be borne in mind in interpreting the results 

which have been given. In any event the restriction does not invalidate the 

conclusions when considered in a comparative sense. 
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(2) The assumption that the drag coefficient is constant demands that the 
"body not he slender and that no shape change occur during entry. We have already 
seen that optimum bodies are, in fact, not slender so that from ttis aspect the 
assumption is valid. However, the process of ablation will tend to promote a 
shape change which, as will be shown later, can and must be minimized so that the 
assumption is therefore valid with this restriction. 

(3) The assumption that ablation is either by sublimation or vaporization of 
the surface is generally an acceptable one for low-temperature ablators but can 
be seriously in error for high-temperature ablators. Quartz, for example, would 
experience considerable ablation in the liquid state unless the heat- transfer 
rates were very high and of short total duration, as during steep trajectories. 
Thus the results given for quartz might be appropriate for the turbulent flow 
cases discussed but surely underestimate energy fraction and mass loss for the 
laminar ones (see ref. 14 for actual performance of quartz at low heat rates). 
Finally, it should be noted that for very high heat -transfer rates, there is con- 
siderable danger that ablation in the solid state may occur as a result of struc- 
tural failure due to excessive thermal stress within the ablator. Stony meteor- 
oids commonly experience structural failure during atmosphere entry (see ref. 20 
and also ref. 21 for a particularly spectacular example) which may result in part 
or in whole from excessive thermal stress. Many of the ablative heat shields 
which have demonstrated excellent performance to date may not fare too well under 
the more severe service we have considered here. 

(4) The assumption is made that the mass lost by ablation during entry is 
small compared to the entry mass. The results show (fig, 13) that Am/m smaller 
than about 0.1 is obtained over a considerable range of entry speeds. For this 
range the assumption is therefore valid. 

(5) The assumption that no cross coupling (ref. 22) of radiative and convec- 
tive heat transfer occurs will certainly be an acceptable one for the results 
which have been given for near-optimum cases since, as has been shown, in these 
cases convective heating dominates. Since the energy loss from the shock layer 
by radiation is small it cannot have an important effect on the convection. At 
conditions away from optimum, coupling may become more important, but such con- 
ditions are not of interest for the present study. 

(6) The assumptions that energy depletion due to radiation is trivial and 
that no reabsorption of radiation occurs within the gas cap are generally accept- 
able for the more interesting case considered (i.e., for near -optimum bodies) 
since the energy fractions are then small. 

( 7 ) The assumption that some limiting Reynolds number exists below which one 
expects laminar flow is open to serious doubt, as has been pointed out earlier. 
This assumption is based on the experiences obtained from tests at far lower 
speeds than are considered in this paper. Clearly, some important research can 
and must be done on the stability of laminar flows at high speeds. A similar 
comment applies to the assumed formulation (eqs. (34) and ( 52 )) for the ratio of 
convective transfer in the presence of ablation to that in its absence. 
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( 8 ) Finally, the assumption that no shape change occurs during entry 
requires considerable review, for it is in contradiction to the fact that the 
heat -transfer rate is not uniform over a conical surface during those portions 
of the trajectory where heating is very important (i.e., in the continuum flow 
regime). In continuum flow, particularly for the laminar case, the convective 
heat -transfer rate varies from large values at the apex of a cone to small values 
at the skirt. Thus the cone tends to he ablated to a round-nosed near-cone with 
increased cone half -angle. If the entry speed is high, the rounded apex, pro- 
moted by the convective heating variation along the cone, becomes flattened by 
the radiative heating contribution at the lower altitudes. The flattened face 
now ablates rapidly because of the near -normal shock conditions at the bow, 
and the mass loss, if the entry speed is high, will be much greater for the 
ablated shape than it would have been if shape change had been prevented. Thus 
such adverse changes in body shape must be prevented if the advantages of conical 
bodies are to be realized. This problem is the subject of the following section. 


The Problem of Shape Change 


To illustrate how the shape of an initially conical body changes with time 
as the result of variations in heat -transfer rate along the surface, a Teflon 
cone was subjected to laminar convective heating in an arc -jet flow of moderate 
enthalpy. (The arc -jet flow characteristics and description of the models used 
in these experiments are given in appendix G. ) Photographs of the body ini- 
tially and at 18 seconds after the establishment of air flow are shown in fig- 
ure 18. The progression of shape change is shown in figure 19* Flattening does 
not occur in this experiment because the enthalpy of the stream is not high 
enough for radiative heating to be significant (the stream enthalpy corresponds 
to a flight velocity less than 3 km/sec). 

If the entry speed of a vehicle is only slightly greater than parabolic 
speed - say 12 km/sec- the penalty for shape changes similar to the initial 
changes shown in figure 19 will, at the worst, be small because, as rounding of 
the apex occurs in an actual flight, the speed diminishes enough that normal 
shock radiative heating, because of its great sensitivity to speed, will not 
cause excessive additional ablation. 

It is anticipated that such will not be the case if the entry speed is some- 
what higher - say Ik km/sec. In this case one might profitably employ a near- 
conical body having a cusped apex, such as shown in figure 20(a). This shape was 
formed by adding to the original cone in the axial direction an amount of Teflon 
very nearly proportional to the inverse square root of the local radius of the 
cone. Figure 20 corresponds to the c = 0.1 case (see appendix G for shape equa- 
tion) . When this cusped shape was subjected to the heating of the arc jet, the 
changes in shape with time were those shown by the photographs and the measured 
ordinates of figures 20 and 21, respectively. One other cusped shape correspond- 
ing to c = 0.3 was also tested. The corresponding changes in shape with time 
are shown by the photographs and measured ordinates of figures 22 and 23 . Com- 
paring the ablation of the cusped bodies with the ablation of the cone, one sees 
that at any given time after the commencement of heating, the nose radius is less 
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the greater the cusp (i.e., the larger the value of c). Thus flight vehicles 
having such cusped shape might he used for somewhat higher entry speeds than 
would he tolerable with a cone. 

The stratagem of cusping cannot he expected to he a satisfactory solution at 
entry speeds well in excess of earth parabolic speed. It was reasoned that for 
such cases one must employ an auxiliary coolant at the apex to prevent drastic 
shape changes. For example, one might cool the cone by transpiration of a gas 
through porous walls at the apex (see, e.g., refs. 23 through 25) or by feeding, 
at an appropriate rate, a fluid or a solid ablator through a hole at the apex. 

To investigate the performance of such an apex cooling system, the Teflon models 
shown in figures 24, 25, and 2 6 were constructed and tested in the arc -jet wind 
tunnel (see appendix G) . The apexes of these models have holes of three differ- 
ent sizes through which Teflon rods of the same diameter as the holes could be 
fed during the tests. Time-sequence photographs of the models during the tests 
are shown in figures 27, 29, and 3T* The model dimension changes with time and 
the feed rates for the rods are given in figures 28, 30, and 32. The remarkable 
fact to note is that this scheme permits the ablation of the conical body to 
occur with but small change in shape of the cone surface. These tests indicate 
that such schemes may be very attractive for application to vehicles at high 
entry speeds. An interesting fact of these particular tests is that the mass of 
rod fed per unit time was essentially independent of rod diameter. This suggests 
that only small diameter rods need be used so that the mass penalty due to rod 
ablation resulting from normal shock radiation at the face of the rod, which 
varies roughly as the cube of the rod diameter, can be kept small. Of course, it 
must be noted that the arc -jet test conditions are considerably different from 
those which would occur in flight and could, therefore, not be representative of 
high entry velocities. The results are nevertheless encouraging that some 
solution to the shape -change problem can be effected. 


CONCLUDING REMARKS 


The analyses employed in this study are clearly of an approximate nature so 
that the results should be regarded as comparative rather than absolute. The 
results do point up future problems of importance for high-speed entries. Some 
salient factors to note are the following: Laminar boundary-layer flows on 

vehicles result in an order of magnitude less mass loss than turbulent flows at 
reasonable reentry Reynolds numbers. Much research on the stability of laminar 
flows at hyperbolic entry speeds must be done in this regard. For optimum coni- 
cal vehicles, convective heat transfer is an order of magnitude more important 
than radiative heat transfer, so that to properly assess heating problems at 
high entry speeds, one must determine convective heating with greater accuracy 
than one need determine the radiative heating. Nonequilibrium radiative heating 
is a much less important factor than equilibrium radiative heating. 

In addition, if the advantages indicated for conical, or near-conical, 
vehicles for atmosphere entry at the very high speeds considered are to be 
realized, changes in heat shield shape due to ablation must be controlled in 


25 



some manner to prevent the serious nose blunting which would normally occur 
during entry. This requirement will demand ingenuity in design. 

Finally, for some entry trajectories of interest (e.g., steep entries with 
turbulent boundary -layer flow) heating rates may greatly exceed those we have 
been accustomed to in the past. Accordingly there is considerable danger that 
many of the ablative materials usually considered to be attractive may fail 
structurally because of excessive thermal stress and so be useless for such 
severe applications. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., July 12 , 1963 


26 



APPENDIX A 


SYMBOLS 


A 

a 

B 


130 ^ 1 ^ 2^3 
C, c 
CD 

Cp 

Ch 

c p 

D 

d 

E 

E 

Ei 

e 

F 

g 

H 

h 

K 

k 


area; without subscript, "base area of cone 

area of shock cone to "base of cone 

ballistic parameter 

constants 

constants 

drag coefficient 

friction coefficient 

over -all heat -transfer coefficient 

specific heat at constant pressure 

drag 

diameter of ablation rod 
kinetic energy 
time rate of radiation 
exponential integral 
Napierian logarithm base 
total frictional force 
constants 

h T 

ratio of total enthalpy, 

h Te 

aerodynamic heat input in kinetic energy units 
enthalpy 

ablative coolant parameter 
thermal conductivity 



I 

M 

m 

Pr 

P 

q. 

R 

Re 

r 

St 

T 

t 

U 

Ui,fi 

u 

V 
v 
X 
x 

Y 

y 

z e 

z 

a 

P 

r 
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model length 
mach number 
mass 

Prandtl number 
pressure 

heating rate per unit surface area 
gas constant 
Reynolds number 
radius 

Stanton number 

temperature 

time 

velocity normal to shock 
a particular velocity, 13,700 ni/sec 
velocity in boundary layer 
vehicle velocity 
volume of shock layer 

distance along axis of revolution measured forward of body base 

distance along the cone surface from apex 

altitude 

distance perpendicular to cone surface 

equilibrium heat-transfer function of sin 

compressibility factor 

function of ablation material 

inverse of scale height 

gamma function 



7 

Am. 
A Q 

5 

i 

n 

e c 

% 

A 

A 

M- 

V 

6 
p 

P 

a 

0 

X 

? 

a 

GJ 


i 


J 1> J 2 


flight -path angle measured from local horizontal 
total ablated mass 

difference between bow shock angle and cone angle 
increment 

energy required to ablate a unit mass 
energy fraction 
half -angle of cone 
shock angle 

function of V E , 0 q, £, and rj 

nondimens ional height within a boundary layer 

absolute viscosity 

kinematic viscosity 

ablation parameter, KV E 2 

air density at any altitude Y 

air density ratio, 

Po 

Ch 

asymptotic value of - — 

C Ho 

function of B 
the product of Bp 
functions of Oq 
functions of £ 
an arbitrary variable 


Powers 


arbitrary 

velocity dependence power in laminar convection 
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ki,k 2 

n 

P 


s 


Id 

cl 

E 

e 

e 2 

fm 

i 

1 

2 1 

lim 

max 

n 

o 

opt 

sp 


velocity dependence power in turbulent convection 
general exponential in power series 
density dependence power in equilibrium radiation 
velocity dependence power in equilibrium radiation 
velocity dependence power in nonequilibrium radiation 
at reference conditions 

Subscripts 

body (or body base as in r^) 
collision limit 
at entrance into atmosphere 
equilibrium radiation 

equilibrium radiation at U < 13,700 m/sec 
equilibrium radiation at U > 13,700 m/sec 
free molecule 
incompressible 
laminar convection 

laminar convection in absence of ablation 
laminar convection for V < 13,000 m/sec 
laminar convection for V > 26,000 m/sec 
limit 
maximum 

nonequilibrium radiation 

sea level 

optimum 

stagnation point 
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turbulent convection 


turbulent convection in absence of ablation 
turbulent convection for V sin < 7^500 m/sec 

turbulent convection for V sin > 7,500 m/sec 

total 

at cone vail 

at edge of boundary layer 
conditions behind the shock 


free stream 


APPENDIX B 


EQUILIBRIUM RADIATIVE HEATING 
Dependence on Velocity and Air Density 


It is By no means obvious that an equation of the form of equation (l 6) can 
he used to represent the dependence on density and velocity of the radiative out- 
put per unit volume of the gas behind a normal shock wave. In particular, if the 
exponents q and p are subject to large and continuous variations and are them- 
selves functions of U and p, the usefulness of this representation would be 
seriously limited. However, it has been previously shown (e.g., in ref. 10) that 
plotting (l/p 1#T ) (dE e /dv) logarithmically against velocity results in a reduction 
of available experimental and theoretical data to a fairly narrow band whose 
center is nearly a straight line for velocities from 5 to 13 km/sec. The theo- 
retical lines in such a presentation show a definite swerving behavior in this 
speed range, indicative of local variations in the value of q (see fig. 7, 
ref. 10), but never deviate too far from a straight line. 

Since we are concerned in this paper with a broader velocity range than that 
considered in reference 10, we plotted the collected theoretical data applicable 
to this broader speed range in figure 33. Here the log 10 (dE e /dv) is plotted 
against log xo U for several altitudes, thus deferring for the moment the con- 
sideration of the effect of altitude. The symbols represent the theoretical data 
and were based on the results of references 5; 26, and 27. It is seen that three 
distinct regions of velocity dependence are suggested by the data, each of which 
can be fitted by a straight line segment. The low velocity range extends from 
about k to 8 km/sec; the intermediate range is from 8 to 13*7 km/sec; and the 
high range extends to at least 30 km/sec. Thus, the two lower speed ranges cor- 
respond to the complete velocity range of reference 10, and the two slopes repre- 
sent, to an acceptable degree, the principal swerve in the theoretical curves 
shown in that reference. In the integrations of total heat input for the present 
paper, the existence of a different velocity dependence in the lowest speed range 
was ignored because of the comparatively small values of radiative heating in 
this speed range; the middle region was assumed to extend to zero velocity. 

In fairing the lines on figure 33* certain restrictions were imposed: 

(l) The line segments in each velocity range were required to have the same slope 
at all altitudes. (2) The intersections of the lines for the upper two velocity 
ranges were required to occur at a constant velocity, independent of altitude. 
These restrictions were made so as to make q independent of p and so that a 
single speed, U ljr2 = 13*700 m/sec, would be the transition speed from one set of 
curves to the other at all altitudes. These restrictions did not seriously 
impair the fit of the lines to the data, the worst error being a factor of 2 and 
occurring at high altitude (and therefore at comparatively low radiative 
intensity) . 

While the theory for the lower -two speed regions is well supported by the 
data of reference 10, the radiative intensities in the highest speed region are 
entirely based on the theory of reference 27 for which there is no experimental 
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verification. This theory treats air at temperatures extending above 200,000° K, 
veil beyond the highest temperatures ( 50 , 000 ° K) considered in the present paper. 
Above 100,000° K, the theory predicts that the radiative energy output is essen- 
tially independent of temperature. Between 20,000° K and 100,000° K, the temper- 
ature dependence may be represented by T 3 ’ 5 . At still lover temperatures, 
according to data given in references 5 and 26 , the dependence is approximated by 
T 1q ' t . Thus the slope change in figure 33 at U l7 2 = 13^700 ni/sec is primarily a 
result of the change of temperature exponent which is predicted to change further 
and go to zero at still higher speeds and temperatures. Although these predic- 
tions should perhaps be viewed with caution pending experimental verification, it 
is noted that analysis of several meteor entries (ref. 28) indicated that the 
magnitude of the radiation predicted by the theory is perhaps correct. 

The altitude or density dependence may now be considered. The radiative 
intensity at a velocity of 13,700 m/sec is plotted on logarithmic coordinates in 
figure 34 for the four altitudes of figure 33 as a function of the density ratio 
p. A straight line is a remarkably good fit to these data, which, together with 
the previous figure, indicates the appropriateness of the form of equation (l6). 
The slope of the line gives p = 1.80. The fact that this plot is made for the 
juncture velocity of the two upper speed regions means that this slope is valid 
for velocities both above and below this velocity. At other velocities, lines 
parallel to that shown are obtained. (if it had not been possible to find a 
common velocity at all altitudes for the intersection of the lines of the two 
families, this would have implied that p is different in the two velocity 
regions.) The evaluation of the constants C ei , Ce 2 , P* q^ an & from fig- 
ures 33 and 34 is straightforward and gives the values recorded in the text as 
equation ( 17 ). 


Calculation of Shock-Wave Angles 


For purposes of calculating the volume of gas radiating to the body and 
defining local flow properties at the boundary-layer edge for the convective 
heating estimates, the following approximate analysis of hypervelocity flow over 
pointed cones was made. The gas was assumed to have the properties of real air 
at equilibrium. The shock layer was assumed to be of uniform density and veloc- 
ity, consistent with the knowledge that the entropy is uniform and the shock 
layer is very thin. Under these assumptions, continuity of flow requires that 

_pV_ 7 + tanAeN 2 ^ = 2 tan A0 + f tan bjf ( B1 

P 2 V 2 \ tan Gq/ tan 0 q \tan 0 q/ 


where p 2 and V 2 are the density and total velocity behind the shock wave. This 
is solved for tan A 0 to obtain 


tan A0 
tan 0 q 


/ — i — 

J 1 - Js OOB e « 


(B2) 
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Given the cone angle 0 q and the free-stream density and velocity, p and V, 
equation (B2) is solved iteratively with the aid of real gas shock-wave tables 
or charts such as those given in reference 29* A first -approximation estimate 
of the standoff angle A0 permits the downstream flow properties to be obtained 
so that pV/p 2 V 2 can be evaluated. Inserted in equation (B2), this leads to a 
new value of AS, and the iteration is continued to convergence (usually one more 
cycle) . 

Standoff angles computed by this technique are shown in figure 35 for veloc- 
ities of 6 to 22 km/sec, cone half -angles of 15°, 30°; and. 55°> and an altitude 
of 50 km. A result of these computations was that the velocity ratio V 2 /V was 
given almost exactly by cos 0^ and within 9 percent by cos 0q, so that equa- 
tion (B2) may be approximated by 


tan A0 
tan 0 q 



(B3) 


or, for p/p £ « 1, A 0 * (l/2)(tan 0 c)(p/p 2 ) radians, a result which is analogous 
to that obtained in studies of shock-wave standoff distance for blunt bodies 
( see, e. g. , ref. 30). 


Since the standoff angle is a function of the density ratio across the bow 
wave, it follows that A 0 will vary in a complicated fashion with free-stream 
velocity and altitude. The most general way to give the results of the present 
calculations is to show the dependence on density ratio, as is done in figure 35, 
where the circular points represent the values obtained at the conditions cited 
in the previous paragraph. For comparison, two points obtained by the exact 
theory of Taylor and Maccoll (see ref. 31) are shown for an ideal gas with a 
ratio of specific heats equal to 1.4 at infinite Mach number (square symbols). 
Since, at the highest speeds considered here, the density ratio at the bow wave 
tends to values between 15 and 18, a working curve was drawn showing approximate 
values of A 0 as a function of cone angle for free-stream velocities greater 
than 20 km/sec, figure 36. For simplicity, this curve was assumed to be univer- 
sally applicable in the equilibrium radiative heating integrations. 

Also shown in figure 3 6 is an interpolated curve from the present calcula- 
tions for a lower speed, 9.2 km/sec, and an altitude of 50 km. This curve is 
compared with points taken from reference 32 for speeds near 9.2 km/sec and a 
free-stream pressure of 0.001 p Q which approximates the pressure at 50 km 
altitude. The method of Taylor and Maccoll for real air at equilibrium was used 
in the calculations of this reference. While agreement is by no means exact, it 
is within 20 percent at the worst point. The reasons for the disagreement are 
not clear. 


A final comment on the behavior of the curves of figure 3 6 in the vicinity 
of 0 q = 0 is in order. Although the lower curve seems to be going into the 
origin, we know that at 0q = 0, A0 = sin -1 (l/M), and the wave angle can never be 
less than the Mach angle (i.e., 00 + A0 > sin _1 (l/M)). Since the Mach angles at 
the speeds considered range from a few tenths of a degree to 1° or 2°, this 
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limitation need not *be considered until < 2°. For cone angles smaller than 

2° (an academic case), the curves must turn up and terminate at A Q = sin -1 (l/M). 


Evaluation of Equilibrium Energy Fraction 

/ . -(B/ 2 )(q, -i)p 

In the integral of equation (22a) the value of e when p is 

of the order of unity or larger is so small that with negligible error we may 
write for p greater than zero 


r 


i -(B/a)(q x -i)p_p-i ~ 

e p ~ 

p=o 


-(B/2)(q 1 - 1 )p.p- 1 
dp = / e p dp = 

o 


r(p) 


B P ( ^ -- 1 
2 


(B4) 


wherein r(p) is the familiar gamma function (see, e.g., ref. 33). Similarly, 
we may write 


'1 -(B/2)(q 1 -i)p p _ x 

e p dp = 


r(p) 


' p i,2 -(B/2)(q 1 -i)p _ p _ 1 


pr dp (B5) 


p=p 


1,2 


B P ( - 1 --— 


wherein 


Pi, 2 = - ln 


2 

B 


-V E sin % 


U 


1,2 


(B6) 


The right-hand integral of equation (B5) may be evaluated by expanding the 
exponential in a power series. Thus, 




s.coi 


-COO i - 

e oo doo - 


to 1 cko = y 

j Q c 1 n * Z_i ni (n + i + 1) 

n=o n=o 


(BT) 


With these formulations, then, the equilibrium energy fraction can be put in 
the form of equation (2k) 


^e z e 


if we put 
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sin 3 ©y(tan 2 ©y - tan 2 © c ) 


sin 2 0 q tan 3 ©q 


(b8) 


Z e = c e x 
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and 


Z e = c ei 
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From the relations of equations (17) and (23), then, 
/V E sin ©y\ 12 * 45 V E sin ©y 

Z e = 0.0701 ( — ^—HL) ; — ^<1.37 
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wherein 


x x = 14.45 2n 


' V E Sin 0 W 

13,700 


x 2 = 4.05 In 


' V E Sln e ¥ 

13,700 


) (Bll) 
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APPENDIX C 


NONEQUILIBRIUM RADIATIVE HEATING 


The treatment of nonequilibrium radiative heating in this paper is based on 
the conceptual model of the nonequilibrium region described in references 34 
and 10, and on the experimental observations of nonequilibrium air radiation of 
reference 10. Although the research on nonequilibrium air radiation is at the 
present writing still very new, it appears to be approaching a satisfactory defi- 
nition of the radiative intensities for speeds normal to the shock wave up to 
13,000 meters/sec. For the higher speeds treated in the present paper, it was 
necessary to extrapolate. 

The data of reference 10 are reproduced in figure 37 on logarithmic 
coordinates as a function of velocity normal to the shock wave. They define a 
line with a slope of approximately 4, and are represented by the equation 


dEn 

da 


1.15X1CT 11 U 4 * os 


The line defined by equation (25a) of the text. 


(Cl) 


— = 0. 74x10 ” 22 U 7 

da 


(C2) 


was adopted on the basis of data available earlier, and is plotted for compari- 
son. For a speed of 11,000 meters/sec, it is seen to indicate about four times 
the radiative intensity given by the latest data. In view of the finding of the 
analysis (see text) that the nonequilibrium radiation makes a relatively small 
contribution to the total heating, no revision to the analysis on the basis of 
equation (Cl) was made. The conservatism implied by use of equation (C2) at the 
higher speeds may be considered desirable in view of the present velocity limit 
of the experiments. 

In spite of the density independence of equations (Cl) and (C2), various 
effects are known to contribute to a cut-off of the nonequilibrium radiation at 
very high altitudes, as discussed in references 3^ en<i 10. The contributing 
effects include the flow energy limitation, truncation of the nonequilibrium 
region by the presence of the body, cooling of the shock layer by a thick viscous 
boundary layer, and collision limiting. For simplicity, these effects are lumped 
in the present analysis into a single effect, referred to as collision limiting, 
and assumed to begin reducing the radiation at altitudes just above 50 km. This 
choice appears reasonable based on information given in the references cited 
above. 


38 



Evaluation of Nonequilibrium Energy Fraction 


The integrals of equation ( 30 ) are 



(C3) 


wherein values of the exponential integral Ei(-w) may be obtained from published 
tables (e.g., ref, 33 ), and 



since is very small. The constant given by equations (26) and ( 31 ) can be 

used to put T] n in the form of equation ( 32 ) 


^n 


= ?n*n 



if we define 


sin 6 % tan 2 0y 
sin 2 tan 2 0 q 
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and 


*n 


3BxlO' T 


Ei(-3B) 


-3Bxio -3 - 

Ei(-3BxlO" 3 ) + 1-6 ^ 

3BX10" 3 J 


(C 6) 


39 


I 


APPENDIX D 


LAMINAR CONVECTIVE HEATING 
Theoretical Method 


The laminar convective heat-transfer coefficients for zero ablation, , 

•were calculated by the method of Kemp, Rose, and Detra (ref. 35) modified as 
described in this appendix. Numerical examples for cone half-angles of 15°; 30°, 
and 55° at velocities of 6, 10, 14, 18, 22 and 2 6 km/sec and at altitudes of 20, 
35, 50, and 75 km were then calculated. With the set of theoretical data so 
obtained, the form of a suitable empirical equation to represent these data was 
determined (eq. (33)) and the empirical constants were evaluated. 

Reference 35 (eq. (21)) gives the approximate relation on which the heat- 
transfer estimates were based which for cones may be written 1 




1 - & 


W 


Pw^w- 


sp 


(Dl) 


where g = h T /h Te , h T = h + (u s /2), h is static enthalpy, g^ = (8g/dA) w , A is 
a transformed coordinate normal to the body surface, e and w indicate conditions 
at the boundary -layer edge and at the wall, respectively, and p, u, and p have 
their usual boundary -layer significance. An approximation suggested in this ref- 
erence is the evaluation of the ratio Pe^e/p^-Mw at the stagnation point for 
subsequent application to all points on the surface. This is said to introduce 
errors no larger than a few percent compared to more exact procedures in several 
instances where they have been compared. Equation (Dl) also assumes that the 
Lewis number is 1.0. 

Since, for a Lewis number of 1.0, 


kw /5h\ 

c Pv W/w 


(D2) 


(where k w and cp^ are the gas conductivity and specific heat at the wall, and 
q is the local heat-transfer rate per unit area), and the relation between y 
and A is as given in reference 35> equation (Dl) may be written for cones as 


q. = 


V3 

2 


0.61*8 







0.438 


sp 


(D3) 


x To avoid conflict with the symbols used in this report, certain of the 
symbols in the equations to follow have been changed from the original form. 
The quantities H and r) of reference 35 are denoted herein as h T and A. 


40 



Integrated over the surface of the cone to obtain dK/dt and divided by 
(l/2)pV 3 A to obtain Ch i 0 , this becomes 


% 


o 


A. 0.6W N/PvUgX ' b/ ^ ( *L) 
43 \ c P w y 


+ (v^/2) - W 

V 2 /2 


P£e 

P W^-W/ sp 


0.438 


(D4) 


where x^ and r^ are the slant length and base radius of the cone and 


\> + (^/ g ) - 

V s / 2 


for high velocities. 

Equation (D4) may be put in terms of the free-stream Reynolds number pVr/p. 
to obtain 


C H l 
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(2/^3)0.648 j p w Mw 
Pr W JpVr b /p ^ p ^ 


p e p € 


^w^w^ sp 


O . 438 


(D5) 


where Pr^ is the Prandtl number at the wall, and u € /V = cos Qq has been used 
(see appendix B) . Properties without subscripts are evaluated in the free stream 
ahead of the bow shock wave. 


Discussion 


In the application of equation (d4), u g was taken to be the value for a 
sharp cone, corresponding to an assumption of small tip blunting. In cases where 
tip blunting occurs to an appreciable extent, u^ becomes lower than it would be 
for a sharp tip, while nothing else in equation (d 4) is affected, so that the 
convective heat transfer is predicted to diminish. In this respect, the results 
of the present analysis are conservatively high. 

The accuracy of the heating estimates obtained on the basis of equation (Dl) 
cannot be assessed at the present time, but they could be in appreciable error. 

It is noted that equation (Dl) is obtained by modifying a stagnation point solu- 
tion on the basis of an empirical observation. Therefore, if it is applied to 
conditions outside the range of its evaluation, as it is here, it may be inaccu- 
rate. Thus, values of ( p G P e /PwMw) S p as low as °-°°35 were encountered in the 

calculations, while reference 35 recommends equation (Dl) for the range of this 
parameter from 0.15 to 0.55* Similarly, values of U£ 2 /h T£ = 2 cos 2 0 q up to 
1.86 were treated, where reference 35 discusses values up to 1-50. Because of 
the uncertainties present, a new theoretical study of real air laminar convection 
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to cones has been started at Ames Research Center, and initial results indicate 
that the present estimates are probably inaccurate at worst by a factor of 2 . 


Computation of Results, and Fitting the Empirical Equation 


The gas properties required in equation (D5) for stagnation temperatures up 
to 47,000° K were taken from references 36 , 37* and 38* Reference 38 gives the 
thermodynamic properties of air for temperatures up to 100,000° K, and refer- 
ence 37 gives the viscosity and other properties up to 15,000° K- Reference 38 
is an extension of reference 37 to 30 * 000 ° K for pure nitrogen, and shows that 
the transport properties do not differ greatly between air and pure nitrogen in 
the temperature range up to 15,000° K. The nitrogen data were therefore used at 
the higher temperatures and extrapolated where necessary. The extrapolation, 
which certainly should not be considered reliable to more than one significant 
figure, was performed on a logarithmic plot of viscosity versus enthalpy. A 
family of curves was plotted with pressure as parameter, so that the shapes of 
all the curves could be used to guide the extrapolation. 

A wall temperature of 3,000° K, representative of the temperatures attained 
by carbon, quartz, and the charring ablators, was selected for evaluation of the 
air properties at the wall. However, it can be shown from equation (D5) that 
Ch 7 is not a sensitive function of the wall temperature. 

L o 

The factor (pe^/pv^Sr^sp 438 e( l’ u - a 'ti° n (D5) is independent of cone angle, 
and is shown plotted in figure 38 as a function of velocity for four altitudes 
from 20 km to 75 km. It decreases continuously with increasing velocity. The 
altitude dependence is small at low velocities, but becomes appreciable at the 
higher velocities. Since an empirical expression in which the density dependence 
is a function of the velocity was undesirably complex for our purposes, and in 
view of the convenience of exponential forms, a single curve independent of the 
density, shown in figure 38 , was taken to represent (Pe^/Pv-Hw-)^! 438 * The equa- 
tion of this line is 


. v 0.438 

/PehA 

VPw^w/ sp 


(3.^xio 6 )v 5/3 


(d 6 ) 


It can be seen from the figure that this expression does not apply at speeds much 
below 10 km/sec, and is very approximate throughout the range. 

The dependence of heat -transfer coefficient on cone angle, according to 

equation (D5), is entirely contained in the radical, J(p^/p) (lo^/fjjcot Qq. The 
factor Pv/p may be expressed 
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which for large values of V sin 6c - U reduces to 


p v V s sin 2 Gq 
P Rz w T w 


(D8) 


and the radical becomes 


/5 5 cot e c - V / 
n/ p m* ° n/; 


2Rz w T w 


sin 20 r 


(D9) 


where J ( ) / 2Rz w Tv is nearly independent of cone angle, speed, and altitude 
and is equal to approximately 0.0016 sec/m for T w = 3000° K. 

Substitution of equations (D6) and (D9) into (D5), along with values for the 
viscosity in the free stream, the wall Prandtl number, and p Q , results in the 
final expression for Ch 7 given in the text as equations (33) and- (38), 
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(DIO) 


Values of C;g 7 given by equation (D5) are plotted against velocity in 
L o 

figures 39 and 40 as points, and are treated as data to be fitted. They show the 
heat -transfer coefficient to be nearly independent of velocity at the lower 
velocities, and then to diminish with further increases in velocity. The curves 
described by equation (DIO) disagree with these data at the low velocities, but 
agree approximately at the higher velocities. Disagreement (in opposite direc- 
tions) also occurs at the highest and lowest altitudes considered (fig. 40). All 
of these trends can be foreseen from figure 38- At speeds below 13 km/sec, a 
reasonable fit is obtained by taking Ch 7 independent of velocity, j = 0 in 

equation (33). The level lines shown in this region correspond to the low-speed 
coefficients of equation (33)* given as equation (37)- ^ may be said that these 

fitted equations describe the calculated data adequately, in view of present 
uncertainties in the method of calculation. 

In applying these equations, a smooth transition was made between the low- 
speed regime with j = 0 and the high-speed regime with j = 1.17* The transi- 
tion curve, described in the text, was a four term series in powers of V, with 
constants selected to match the levels and slopes of the two heating equations 
at 13,000 m/sec and 26,000 m/sec, respectively. 



Evaluation of Laminar -Convective Energy Fraction 


To solve equation (41+) we make the following assumptions: 

(a) Because when p = 1, the exponential 

e - B P (Dll) 

is very nearly zero for vehicles of interest to this paper, it is satisfactory to 
substitute unity for 


1 + | z e' B5 (D12) 

(h) Because in the range 0 < p < p „ for vehicles of interest to this paper 

in. — B o 

v = Vg, then for the denominator of the first integral the expression 1 + M 

may be replaced with 1 + % ^ . 

With these simplifications equation (44) becomes, if we let 
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The first and third integrals can be evaluated numerically; the second is a 
gamma function; the fourth can be evaluated as a series after expanding the expo- 
nential as a series; and, although the fifth can be integrated directly, it is 
preferable to express it in series form also. The solution may be written as 
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Since (see ref. 33) r(3/2) = 0.8862, then equation (D15) may he -written 
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For vehicles of interest in this report the finao. term shown above in this 
equation is always negligibly small but the preceding term can be important for 
"small light" entry bodies (very large The value of <6 Z is only impor- 

tant for "large heavy" entry bodies (very small B). Values of 0i 1 and 0i £ , 
and and $i £ (corresponding to j x and j 2 of eqs. (37) (38)) are given 

in tables IV and III, respectively. 
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Equation (Dl8) can be -written in the form of equations (45) and (46) if we 


set 


? 7 , = ii = 1.76x10 

1 2 sin 2 0 q 


_ 4 s/sin 200 


sin 2 00 
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_ 4 s/sin 200 


2 sin 2 00 sin 2 00 

Values for these functions are given in table I. 
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APPENDIX E 

TURBULENT CONVECTIVE HEAT TRANSFER 
Method of Estimation 


Turbulent heat -transfer rates for pointed cones at flight speeds up to 
30,000 m/sec were estimated by use of Colburn's modified Reynolds analogy, ref- 
erence 39> and the reference enthalpy method. The latter is a generalization to 
the case of a dissociating gas at high temperature of the reference temperature 
or T' method described in reference 40. The T T method has been compared with 
measurements of skin friction and heat transfer for turbulent boundary layers in 
subsonic, supersonic, and hypersonic flow of air in or near the ideal gas temper- 
ature range, and is found to represent existing experiments unusually well (see, 
e.g., refs. 4l and 42). Of course, no experience exists for the much higher 
speeds considered here, and, in light of the purely empirical nature of the 
method, the degree of accuracy of these estimates cannot be presently determined. 
However, the excellent experience at supersonic Mach numbers below 10 encourages 
us to expect realistic values of heat -transfer prediction in the higher speed 
range. 

The T' equation of reference 40 is 


|^ = 1 + 0.035 M € s + O .45 (h - 1 ) (El) 

which may, for the ideal gas considered in that reference, be written in the form 

-1 + 0.175 (4$ - 1)+ 0.1*5 (g - 1) (E2) 

where T S p is the stagnation point temperature. These equations define the 
temperature T* within the boundary layer at which to evaluate the density and 
viscosity so that the incompressible skin friction equation is applicable to the 
compressible flow boundary layer considered. The generalization to high tempera- 
ture flow of equation (E2) is assumed to be 

ET 1 + °- 175 (iff - x ) + °-‘* 5 (s' 1 ) <E3) 

which reduces to equation (E2) for a gas with constant specific heat. Here, 
(h S p/h G ) - 1 = (l/2) (u € 2 /h e ) . 1 Given h' as defined by equation (E3) and the 

Although the details will not be given here, it can be shown starting from 
equation (E 3 ) that the reference enthalpy is approximately a fixed fraction of 
the stream total enthalpy for any cone angle, the fraction being a direct 
function of the cone angle. 
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local static pressure on the cone, the state of the gas at the reference enthalpy 
condition is fully specified and the reference density p' may he obtained from 
reference 36 and the reference viscosity p' may he obtained from reference 37 
or 38. The Reynolds number Re' given by p'u e x/p' is then inserted in the 
incompressible turbulent skin -friction formula to obtain Cp', which is related 
to the average skin friction over the surface Cp by 

Cp | p e u £ £ = Cp' | p'u e 2 


or 


CpP € - Cjp 1 p ! 


The variation with Reynolds number of average skin friction coefficient on 
a flat plate in incompressible flow was taken to be that given by the K&rm&n- 
Schoenherr formula 


~= = log 1Q ReCp (El) 

which can be very closely approximated for Reynolds numbers from 3X10 6 to 10 9 by 
the power law formula 



0.0317 

O. 148 
Re 


(E5) 


The turbulent skin friction values for a flat plate were modified by the factor 
1.047, derived in appendix C of reference 39 > to make them applicable to cones. 

The Colburn modification of the Reynolds analogy 


St 


Cp 

ss 7 73 


(e6) 


was then applied to obtain the dimensionless heat -transfer parameter, the Stanton 
number, defined by 


St (ET) 

Pe u e(hsp - hw)A w 

where A v is the wetted area of the cone. From dH/dt, Cg. may be obtained 

^o 

from its definition, equation (4) of the text. 

The above steps describe the procedure used, 
mating these steps may be written 


A summary equation approxi- 



(E8) 



0.0166 

/ . , . . 0.148 

(p'ueXb/n' ) 


p . cot e c 
P Pr w 2/ 3 


in which u e /V = cos Gq and (h S p - h^) /V 2 ~ l/2 have been applied along with the 
approximation given hy equation (E5) for the KdrmAn -S cho enherr formula. 


Expression in Terms of Free-Stream Variables 


Equations (e 8) and (E3) describe, in reasonably simple terms, the heat- 
transfer coefficient for the turbulent boundary -layer case, but they are not 
suitable for incorporation in the optimization analysis of this paper because 
they are expressed in terms of p T and \±' , the variation of which with speed, 
cone angle, and free-stream density is not immediately evident. The point of 
departure for finding ^E±, 0 as a function of V, 0, and p was not equation (E8), 

but was based on the observation that for any altitude and base radius, the 
Stanton number was a nearly linear function of 0 q, almost independent of veloc- 
ity. This is shown in figure 41 for an altitude of 35 km* The equation for 
Stanton number obtained from plots of this kind for altitudes of 20, 35, and 
50 km is 


St 


2.47 9 0 

nn 5 ™ 0 . 148 - 0.165 

10 r t) P 


(E9) 


where 0 q is in degrees and r^ in meters. The fact that the exponent of p 
is here found empirically to differ a little from 0.148 reflects the dependence 
of the Stanton number on p'/p and M-’/m* which vary with altitude for any given 
cone angle. 

From equations (E7), (E9), and (4) in the text, 

2.47 0 C P6 u e^sp “ ^w)^w 

H to = 10 5 Pb° ■ ■ 14s p° ‘ 165 ( 1/2 ) pV 3 A 


2.47 e C p e 

" 10 5 1^°- 148,0. 165 p co1: y C 


(E10) 


where u € /V = cos 0g, h S p - h w » V 2 /2, and A w /A = l/sin 6q have "been applied. 

In addition, the density ratio across the shock wave, p e /p, may he correlated in 
terms of the velocity component normal to the shock wave, V sin 0^ = U, as shown 
in figure 42. Two regions exist. For U < 7500 m/sec, the density ratio 
increases with increasing speed approximately as 


= ~ 5?039 ( V Sln ®W) 


O . 66 


(Ell) 
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For U > 7500 m/sec, the density ratio becomes a slowly decreasing function of 
increasing velocity. 


p e 54.5 

p pO • 039 


(V sin %) 


(E12) 


The accuracy of these expressions in relation to calculations for normal shock 
waves in equilibrium air may be examined in figure 42. Inserting these relations 
into equation (E10) leads to the expressions for given as equations (50) 

and (51) of the text. 0 

The values of computed from the working equations (50) and (51) are 

compared in figure 43 for a cone angle of 30° and a 1-meter base radius with the 
values computed by the complete estimation procedure based on equations (E3) 
and (e 6) . The representation given by the working equations is generally ade- 
quate although, for this cone angle, uniformly a little high. The values of con- 
vective heat -transfer coefficient estimated for a laminar boundary layer on a 
30° cone are reproduced on this figure for comparison and are computed to be 
smaller than those for a turbulent boundary layer by factors ranging from 3 to 
40, depending on the speed and altitude. 


Evaluation of Turbulent Convective Energy Fraction 


To solve equations (53 a ) and (53b) we make the assumption that when p = 1 
the exponential e _ ®P is so small that we may substitute unity for (l + . 

With this simplification equations (53 a ) and (53b) "become 

for 


y < 7500 

sin % 


^t 


^0.204 

l"c+ 

tl 1 

f Qq cos Qq \ 

y. 0.148 

r b 

V % 

LV e 1 

^2 sin* 1 9\fj sin 3 QqJ 






dX 


Ct W 

' 0q cos 

0 c \ 

r axi 


v 2 sin^ 1 Q\\ 

sin 3 00^ 

J B x 0 - 204 J 


and for 

v > 7500 

sin 


(E13) 
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•gO . 204 
= „ O. 148 


Ct x / © C cos e c 


V E kl k 2 sin kl % sin 3 x= 2 ln( 


X= 2 Znl^| i^) Vl ^ te 


1 ' at , „ e 
^( +CT t 


ax 

^0.204 


ctfe ^ &c cos e C 


V E k2 ^2 sin^ % sin 3 ©0 7 w ° 
c t x f ®C cos 9 C 


. 


:)/ 


/ V E sin %\ 

\ 7500 / ✓ 1 


1 + I t e 


-X 


+ °t 


^ ax 

v 0 . 204 


V E 1 ^2 sin 1 ©VJ sin 3 ©c 7 


-(w^MX .1 


<3X 


,0.204 


(ElU) 


where 


X = Bp 

In equation (El4), the first two integral terms, in addition to Being functions 
of at and £t> are functions of Oq and in a way which is not easily sepa- 
rable. For the purposes of this paper it was deemed advisable to leave the sum- 
mation in this inseparable form so that a two -entry table is required for tabula- 
tion, and each table, in turn, must be made up for a particular value of at and 
Accordingly, we may put in the form of equation (55) 


^t = 


B' 


0.204 
O. 148 


A t - H*t 


Ve\ 

lo v 


0.66-1 


if we set 


*t = I ^xio" 43 " 1 


©Q COS ©0 


sin kl % sin 3 © E - 


= 0. 202X10 ' 3 


©0 sin 0, 66 0 y cos © 0 ' 
sin 3 00 


(E15) 


H = 


"°° e -°- 67y ax 


X 

while for 


0 . 204 


v < 1^22 
sin ©^ 


(El 6 ) 
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and for 


V> 


7500 

sin % 


At = 




2V E kl Vsn-n 1 ^ 1 


0 Q COS 00 


1 - 




-sin 1 0w sins 9? x=sln ( - E sln + | t e 

\ 7500 / 

c t 2 f 9 C cos e c 


-X 
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2Vg k2 Nsin^ 2 sin 3 o 
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APPENDIX F 


RELATION OF THE LIMITING REYNOLDS NUMBER TO THE BALLISTIC 
COEFFICIENT, CONE ANGLE, AND ENTRY VELOCITY 

Since the relationship Between flight velocity and air density is a function 
of the entry velocity and Ballistic coefficient B (eq. ( 8 )), it is evident that 
the Reynolds number Based on free -stream properties will Be governed By these 
factors and the Body size. 


PVT, _ -(B/*)p 

00 \i q 


(Fl) 


For an isothermal atmosphere in which the free-stream viscosity is constant, and 
for cases where is not importantly diminished By aBlation, the maximum value 

of Reoo on the trajectory can Be found By setting the derivative dRe^/dp equal 
to zero to obtain 


_ _ 2 _ Pp V E r l3 
e °°max Be n 


(F2) 


where e is the Napierian Base. Equation (F2) is plotted in figure 44 for an 
entry velocity of 30 km/sec and r^ = 1 meter, for which case values of Reoomax 
in excess of 10 million occur at B less than about 265- It is of interest that 
the peak Reynolds number according to equation (Fl) occurs at a velocity 
V = Vg/e, and at a density p = 2/B. The Reynolds number profile defined By 
equation (Fl) is plotted in figure 45 as a function of the velocity ratio V/Vg. 

It remains to relate the free-stream Reynolds number defined By equa- 
tion (Fl) to the Reynolds number Based on Boundary-layer edge properties and 
slant length, Re e = p £ u e x/q £ . All flow-field properties needed for this purpose 
are obtainable from the calculations of the conical flow field described in 
appendix B, with the use of air properties tabulated in references 36 , 37, 
and 38 . Ratios of local to free-stream Reynolds number obtained from these cal- 
culations are shown in figure 46. 

Since the present purpose was to define limiting values of B Below which 
the local Reynolds number will exceed 5, 10, or 20 million, the following pro- 
cedure was applied. Assume for the moment that the Reynolds number ratios of 
figure 46 are functions of cone angle only, independent of speed and altitude. 
Then the maximum of the local Reynolds number will occur simultaneously with the 
maximum of the free-stream Reynolds number, and 


Ree max “ Re max - Re oo m ax 


Ree _ 2 ^o^E r b R e e 
Reoo Be M' Re^ 


(F3) 
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If Re e max is not to exceed the limiting Reynolds number Re lim of the text, 
then the former is replaced in equation (F3) hy the latter to give the value 


2 p o V E r h Reg 
eRe lim ^ Re °° 


(F4) 


corresponding to the specified Reynolds number limitation. Lower values of B 
will result in local Reynolds numbers exceeding Req im . 

As may he seen from figure 46, the assumption that Re e /'Re O0 is a constant 
for any cone angle is reasonable for cone angles up to 40° at speeds up to 
25 km/sec, and to higher speeds for smaller cone angles. Obviously, large depar- 
tures from constant Re 6 /Re^o occur on the 40° and 55° cones at speeds above 
20 to 25 km/sec. The error resulting from the assumption of constant Re e /'Re Q0 
in this latter range is not important for present purposes, however, because the 
optimum cone angle in the very high entry speed range is less than 33 0 (see 
fig. 10). Hence, the value of B calculated from equation (F4) by use of this 
assumption will be in serious error only for conditions well away from the opti- 
mum, and the effect on conditions near the optimum will be small. It is of 
interest to note that the effect of the nonconstancy of Re e /Re b0 on the dashed 
curves of figure 9 is to cause them to rise more rapidly than shown for 0q 
greater than the optimum. The significance of this is that the turbulent flow 
boundary is raised to higher values of B for the large angled cones. Since 
these cone angles would presumably not be used, this may be of academic interest 
only. It is clear that by use of equation (Fl) and figure 46, values of maximum 
local Reynolds number free of the above assumption can be obtained if desired. 
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APPENDIX G 


ARC -JET TESTS OF TEFLON MODELS 


Tests to examine qualitatively the extent of the change in dimensions of 
Teflon models due to laminar convective heating were made in an electrically 
heated arc -jet stream. The Mach number of the stream was 3.3 and the enthalpy 
level, 3500 Btu/lb (8.1xl0 6 m 2 /sec 2 ). Total pressure was set at one atmosphere. 
The Reynolds number based on model base diameter was 2140 and, accordingly, the 
flow was laminar. The jet diameter was 2.75 inches (6.8 cm). Stream surveys 
indicated the enthalpy and stagnation pressure were constant within 5 percent 
over a 2-inch core (5 cm). This assured uniform flow conditions over the conical 
test bodies which were 1.5 inches (3-8 cm) in diameter. 

The test conditions in this stream are far less intense, of course, than a 
typical flight vehicle would experience at hyperbolic entry speeds (e.g., in 
flight with a ballistic parameter of 200, heat rates would be as much as 50 times 
those in these tests) and the jet enthalpy is so low that radiative heating of 
these models is trivial. However, the model test results should be indicative of 
the manner in which convective heating promotes the initial blunting of the apex 
in flight. 

For the flow conditions listed above for the present ablation tests, the 
following table lists the magnitude of the stagnation -point heat -transfer rates 
at 1 second after flow was started. 


Model designation 

Btu/sec 

ft* 

q k S 

sec 2 

Passive type; c = 0 

930 

1.05X10 7 

c = . 1 

1440 

1.63XIO 7 

c = .2 

1660 

i.88xio 7 

c = .3 

1850 

2.ioxio 7 

Extrusion type ; l/l6 -diameter rod 

930 

i.05xio 7 

l/8-diameter rod 

930 

i.05xio 7 

l/4-diameter rod 

930 

1.05X10 7 


Figure 47 shows the principal dimensions of the passive type models tested. 
Basically, all of the passive type models are members of the same cusped-cone 
family. The degree of sharpness of the nose of the cusped cone is determined by 
the parameter c (see equation, fig. 47)- For practical reasons of strength, all 
models were terminated at a tip diameter of l/64 inch. The notch on the cylin- 
drical afterbody served as a measuring reference station. 

The other type of model tested is shown in figure 48. This model consists 
of a 45° cone with a cylindrical core of various diameters that could be extruded 
at constant predetermined rates. Extrusion rates were set by adjusting a D.C. 
motor- driven drive screw prior to each run. An "optimum" extrusion rate was 
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found "by trial. The optimum rate was defined such that the core tip just kept 
pace -with the ablation that occurred on the conical rays. For this case, the 
profile shape remains essentially constant with time throughout the run. 

The profile shape change with time was found for all models by measurements 
from shadowgraph pictures taken during the run. 
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TABLE I.- FUNCTIONS OF 0 C 
[Code: 0.76-n = 0. 76x10 “ n ; 0.76+n = 0.76xl0 n ] 


ec 

sin 


*e 



**2 


15 

0.267 

0.13 

O.76-1 

0.58-2 

1.86-3 

2.18-3 

0.67-1 

20 

.354 

• 23 

.83-1 

I.8I-2 

1.21-3 

1.42-3 

.45-1 

25 

.438 

.36 

• 90-1 

.43-1 

.85-3 

1.01-3 

.33-1 

30 

.519 

.50 

.98-1 

.81-1 

.65-3 

.77-3 

.26-1 

35 

.595 

. 66 

1.07-1 

1.50-1 

.52-3 

. 61-3 

2.03-2 

40 

.667 

.82 

I.I8-1 

2.41-1 

.42-3 

.50 -3 

1.65-2 

45 

•734 

1.00 

1.31-1 

•37 

.35-3 

. 4-1-3 

1.35-2 

50 

•795 

1.17 

1.47-1 

• 52 

.29-3 

.35-3 

1 . 11-2 

55 

.848 

1.3^ 

I.63-1 

•70 

. 25-3 

.30-3 

.91-2 


TABLE II.- FUNCTIONS OF VE sin 0 W 
[Code: 0.76-n = 0. 7^X10 -n j 0.76+n = 0 . r j6xlO n ; V E in meters/sec] 


Ve sin 

z e 

V E sin e w 

Ze 

6,000 

2.40-4 

18,500 

. 68+2 

6,500 

.66-3 

19,000 

. 76+2 

7,000 

1.66-3 

19,500 

. 84+2 

7,500 

.39-2 

20,000 

. 94+2 

8,000 

.88-2 

20,500 

1.02+2 

8,500 

1.86-2 

21,000 

1.11+2 

9,000 

.38-1 

21,500 

1.20+2 

9,500 

.74-1 

22,000 

1.29+2 

10,000 

1.40-1 

22,500 

1.39+2 

10 , 500 

2.58-1 

23,000 

1. 48+2 

11,000 

.46 

23,500 

1. 58+2 

11,500 

.80 

24,000 

1.68+2 

12,000 

1.36 

24,500 

1.77+2 

12,500 

2.26 

25,000 

1.87+2 

13,000 

.37+1 

25,500 

1.97+2 

13,500 

.591-1 

26,000 

2.07+2 

14,000 

.92+1 

26,500 

2 . 16+2 

14, 500 

1-33+1 

27,000 

2.25+2 

15,000 

l. 86+1 

27,500 

2.34+2 

15,500 

2. 44+i 

28,000 

2.43+2 

16,000 

.31+2 

28,500 

2. 50+2 

16,500 

•37+2 

29,000 

2 . 5 6+2 

17,000 

.46+2 

29,500 

2.64+2 

17,500 

18,000 

.52+2 

.60+2 

30,000 

2. 691-2 
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TABLE III. - FUNCTIONS OF B 
[Code: 0.76-n = 0 . 76XIO _n ; 0.76+n = 0 . 76 xl 0 n ] 


B 


$ Zi 

$ Z 2 

H 

l B 


B 

4 >n 

0 

0 

1.77 

2.75 

1.56 

25 

2 . 29-5 

160 

. 66-4 

1 

1 . 87-6 

2.79-1 

1.00 

.66 

30 

2 . 59-5 

170 

.68-4 

2 

• 33 -s 

.81-1 

•54 

■31 

35 

2.86-5 

180 

.70-4 

3 

. 46 - 5 

2 . 54-2 

• 32 

1.49-1 

40 

• 31-4 

190 

.72-4 

4 

.58-5 

.83-2 

1.88-1 

.73-1 

45 

• 34-4 

200 

.73-4 

5 

.69-5 

2.78-3 

1.15-1 

.36-1 

50 

.36-4 

250 

•79-4 

6 

.80-5 

. 94-3 

.71-1 

1.78-2 

55 

.38-4 

300 

.83-4 

7 

.90-5 

• 32-3 

. 44-1 

.89-2 

60 

. 40-4 

350 

.86-4 

8 

• 99-5 

1 . 12-4 

2 . 74-2 

. 44-2 

65 

. 42-4 

400 

.89 “4 

9 

1.09-5 

• 39-4 

1.73-2 

2.23 -3 

70 

. 44-4 

450 

.91-4 

10 

1.18-5 

1 - 37-5 

1.09-2 

1 . 12-3 

75 

. 46-4 

500 

• 93-4 

11 

1.27-5 

. 48-5 

.69-2 

.56-3 

80 

.47-4 

600 

• 95-4 

12 

1.36-5 

I.71-6 

• 44-2 

2.83-4 

85 

.49-4 

700 

! - 97-4 

13 

1 . 44-5 

.61-6 

2.81-3 

1 - 43-4 

90 

.50-4 

800 

. 98-4 

14 

1.52-5 

2 . 15-7 

1.80-3 

• 72-4 

95 

.52-4 

900 

.98-4 

15 

1.60-5 


1.15-3 

• 36-4 

100 

• 53-4 

1,000 

.99-4 

16 

1.68-5 


.74-3 

1 . 84-5 

110 

. 56-4 

2,000 

1 . 00-4 

17 

1.72-5 


.47-3 

• 93-5 

120 

.58-4 

3,000 

1.00 -4 

18 

1.82-5 


.31-3 

.47-5 

130 

.60 -4 



19 

1.89-5 


1 . 97-4 

2 . 39-6 

140 

.62-4 



20 

1.96-5 


1.27-4 

1 . 21-0 

150 

. 64-4 



TABLE IV.- FUNCTIONS OF 





[Code: 0.76-n = 0.76xl0” n ; 0.76+n = 0.76xl0 n ] 


h 

a lx 

a l 2 


n z 2 

0 

1-77 

2-75 

90 

2 . 42-1 

5 

.46 

1.04 

100 

2.29-1 

10 

2.82-1 

•75 

no 

2.18-1 

15 

2.06-1 

.61 

120 

2.09-1 

20 

1.64-1 

•53 

130 

2.00-1 

25 

1 - 37-1 

•47 

i 4 o 

1 . 93-1 

30 

I.18-1 

•43 

150 

1.86-1 

35 

l. 04 -i 

•39 

200 

1.61-1 

40 

• 93 -i 

•37 

250 

1 1.43-1 

45 

. 84 -i 

• 35 

300 

1.31-1 

50 

.76-1 

• 33 

350 

1.21-1 

55 

.71-1 

• 31 

400 

1.13-1 

60 

.65-1 

2.98-1 



65 

.61-1 

2.86-1 



70 

■ 57-1 

2.76-1 



75 

• 54-1 

2.66-1 * 



80 

.51-1 

2 . 57-1 1 





TABLE V.- VALUES OF At 
[Code: 0.76-n = 0.76xl0 _n ; 0.76+n = 0.76xl0 n ] 


(a) 



m/ sec 

15° 

8,000 

0.53-1 

10,000 

•53-i 

12,000 

•52-1 

14,000 

.51-1 

16,000 

50-i 

18,000 

■ 50-1 

20,000 

• 49-1 

22,000 

.49-1 

24,000 

• 49-1 

26,000 

.48-1 

28,000 

.48-i 

30,000 

.48-i 


(a) Subliming Teflon ablator, Kt = 4x10 ~ 8 sec z /?(£ , a t = 0.1 

0c 


20° 

25° 

O 

O 

CO 

35° 

40° 

O 

Lf\ 

-3“ 

50° 

55° 

36-1 

2.61-2 

2.01-2 

1 . 60 -2 

1.30-2 

1.07-2 

0 . 88 -2 

0.72-2 

35-i 

2.60-2 

2.00-2 

1-59-2 

1.29-2 

1.06-2 

.87-2 

.71-2 

35-1 

2.56-2 

1.98-2 

1-57 -2 

1.27-2 

1.04-2 

.85-2 

.69-2 

35-1 

2.53-2 

1.95-2 

1.54-2 

1.24-2 

1.01-2 

.82-2 

.68-2 

34-1 

2.49-2 

1.91-2 

1.51-2 

1.21-2 

. 98 -2 

.79-2 

.64-2 

34-1 

2.46-2 

1.87-2 

1 . 47 -2 

1 . 17 -2 

.94-2 

• 77-2 

.62-2 

33 - 1 

2.42-2 

1.83-2 

1.43-2 

1.14-2 

.92-2 

.75-2 

.61-2 

33 - 1 

2.38-2 

1.80-2 

1.40-2 

1.11-2 

.90-2 

•73-2 

• 59-2 

33 - 1 

2.34-2 

1.76-2 

1-37-2 

1.09-2 

.88-2 

.71-2 

.58-2 

32-1 

2. 31-2 

1-73-2 

1.34-2 

1.07 -2 

.86-2 

.70-2 

.56-2 

32-1 

2.28-2 

1.70-2 

1.32-2 

1.05-2 

.84-2 

.68-2 

• 55-2 

32-1 

2.25-2 

1.68-2 

1.30-2 

1.03-2 

.83-2 

.67-2 

• 54-2 


(b) Vaporizing quartz ablator, Kj- = 0.5X10 -8 sec 2 /m 2 , ct = 0*1 


.82-1 

.90-1 

.96-1 

.99-1 

1 . 02-1 

1.03- 1 

1.04- 1 

1.05- 1 
1.05-1 
1.05-1 
1.04-1 
1.04-1 


2 . 50-2 
2.73-2 
2.9O-2 
3.OO-2 
3.04-2 
3-03- 2 
3.00 -2 
2.97-2 
2.92-2 
2.88-2 
2.83-2 
2.78-2 


2.03-2 

2 . 22-2 

2 . 35-2 

2.41- 2 

2.42- 2 
2.41-2 
2.38-2 
2-35-2 
2.31-2 
2.27-2 
2.24-2 
2 . 20-2 


1.66-2 

1.82- 2 
1.91-2 
1.95-2 
1-95-2 

1.94-2 

1.91-2 

1.89-2 

1.86-2 

1.83- 2 
1.79-2 
1.76-2 


1-37-2 

1.49-2 
1.56-2 
1.58-2 
1.58-2 
1.56-2 
1.54-2 
1.52-2 
1.50 -2 
1 . 47 -2 
1.45-2 
1.42-2 


1 . 12-2 

1 . 22-2 

1.26-2 

1.28-2 

1.27-2 

1.26-2 

1 . 25-2 

1.23-2 

1 . 21-2 

1.19-2 

1.17-2 

1.15-2 
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Figure 3 .- Variation of the laminar energy fraction with entry speed for a 
30° half -angle cone with a ballistic parameter of 200 and a Teflon 
ablator. 



Figure 4.- Variation of the laminar energy fraction with entry speed for cones 
having a ballistic parameter of 200, Teflon ablators, and various cone 
half -angles. 
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Figure 9 »- Example of a plot for obtaining optimum laminar energy fractions 
and allied factors for three limit Reynolds numbers. 
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Figure 11.- Comparison of optimum energy fractions of Teflon and vaporizing 
quartz for a limit Reynolds number of 10 7 . 



Figure 12.- Generalized mass-loss factor as a function of entry speed. 






Figure 15* - Trajectory angle requirements for maintenance of laminar flow for 
an entry body having typical space -probe characteristics. 



Figure 1 6 .- Variation of minimum energy fraction with ballistic parameter for 

several entry speeds. 
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Figure 17*- Variation of cone half -angle corresponding to minimum energy 
fraction with "ballistic parameter for several entry speeds. 





A-30693 


Figure 18.- Conical "body before and after exposure to arc-jet stream; 
h T = 35 OO Btu/lb; p S p = 0.114 atm; M = 3*3J exposure time = 18 sec. 
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A-30694 


Figure 20.- Cusped-cone shape (c = 0.1) before and after exposure to arc-jet 
stream; h T = 3500 Btu/lb; p S p = 0.114 atm; M = 3*3* exposure time = 18 sec. 
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A-30842 


Figure 22.- Cusped-cone shape (c = 0.3) "before and after exposure to arc-jet 
stream; h T = 3500 Btu/lb; p S p = 0.114 atm; M = 3-35 exposure time = 18 sec. 
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Figure 24.- Conical "body shape (45°) with a l/l6-inch- diameter 

extruded rod. 



A-30818 

Figure 25-- Conical body shape (45°) with a l/8- inch- diameter 

extruded rod. 
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Figure 26.- Conical "body shape (k5°) ‘with a l/4- inch- diameter 

extruded rod. 
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A-30695. 1 


Figure 27.- Conical body with a l/l6-inch-diameter rod. extruded during exposure 
to arc -jet stream; h T = 3500 Btu/lb; p S p = 0.114 atm; M = 3«3j exposure 
time = 18 sec; feed rate = 0.057 in/sec. 
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Figure 28.- Profile shape change during ablation for a ^5° cone with a l/l6-inch 
diameter extruded core; base radius = 0.75 inch; time interval = 6 sec; 
feed rate = 0.057 in/sec. 





A-30840 


Figure 29.- Conical "body with a l/8 -inch-diameter rod extruded during exposure 
to arc-jet stream; h T = 3500 Btu/lb; p sp = 0.114 atm; M = 3-3; exposure 
time = 18 sec; feed rate = 0.034 in/sec. 
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Figure 30.- Profile shape change during ablation for a 45 ° cone with a l/8-ineh 
diameter extruded core; base radius = 0-75 inch; time interval = 6 sec; 
feed rate = 0.034 in/sec. 





A-30839 


Figure 31.- Conical body with a l/4-inch-diameter rod extruded during exposure 
to arc -jet stream; h T = 3500 Btu/lh; p sp = 0.114 atm; M = 3.3; exposure 
time = 18 sec; feed rate = 0.024 in/sec. 
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Figure 32 .- Profile shape change during ablation for a 45° cone with a l/4-inch 
diameter extruded core; base radius = 0-75 inch; time interval = 6 sec; 
feed rate = 0.024 in/sec. 







U 7 meters/sec 

Figure 3 ?.- Nonequilibrium radiation intensity as a function of Telocity 

normal to the shock wave. 



m/sec 

Figure 38.- Variation of the density-Tiscosity factor as a function of 

altitude and speed. 






Figure 39 .- Variation of heat -transfer coefficient with speed at 50 km altitude 
for several cone angles and a "base radius of 1 meter. 




Figure 40.- Variation of heat-transfer coefficient with speed for 30° cone angle 
at three altitudes for a base radius of 1 meter. 
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Model designation rt> cr b 2 

c - 0 0.75in (1.9cm) Oin (Omm) 0.75in (190cm) 

£ = 0.1 0.75 in (1.9cm) 0.075in (1.9 mm) 1.08 in (2.74cm) 

£=0.2 0.75in (1.9cm) 0. 15 i n (3.8mm) 1. 45 in (3.67cm) 

£=0.3 0.75 in (1.9cm) 0.225 in (5.7 mm) 1-81 in (4.58cm) 

Figure 47 *- Cusped-cone models. 



Model Core 

designation diameter, d 

No. I 0.25in (6.35mm) 

No. 2 O.I25in (3.17 mm) 

No. 3 0.0625in (1.58 mm) 

Figure 48 . - Extruding apex models . 
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